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Abstract 



Phase and modulus of an energy- and pressure- free, composite and adjoint field 
in an SU(2) Yang-Mills theory are computed. This field is generated by trivial 
holonomy calorons of topological charge one. It possesses nontrivial S\ -winding 
on the group manifold. The two-loop contribution to the thermodynamical pres- 
sure of an SU(2) Yang-Mills theory in the electric (deconfining) phase is com- 
puted in the real time formalism of finite temperature field theory. The result 
supports the picture of only very weakly interacting quasiparticles. 



Zusammenfassung 



Phase und Betrag eines energie- und druckfreien, zusammengesetzten und ad- 
jungierten Feldes in SU(2) Yang-Mills Theorie werden berechnet. Dieses Feld 
wird von Kaloronen trivialer Holonomie und topologischer Ladung eins erzeugt. 
Es besitzt nichttriviales Verhalten auf der Gruppenmannigfaltigkeit. Der Zwei- 
Schleifen Beitrag zum thermodynamischen Druck einer SU(2) Yang-Mills The- 
orie in der elektrischen (deconfinierten) Phase wird im Real-Zeit-Formalismus 
berechnet. Das Ergebnis stiitzt das Bild nur sehr schwach wechselwirkender 
Quasiteilchen. 
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Chapter 1 
Introduction 



Gauge theories such as the Standard Model of particle physics have been investigated mainly 
in the framework of perturbation theory. Due to the tremendous complexity of (especially 
nonabelian) gauge theories, this is the only feasible approach for problems concerning single 
fundamental processes: It is not possible to solve even the simplest QED processes exactly. 
Concerning the macroscopic properties (e. g. thermodynamics) of a system, chances for a 
nonperturbative treatment in form of an effective theory are much better. Superconductivity 
in metals, for example, was described phenomenologically by the Landau- Ginzburg theory PU 
12] . Subsequently, a quantum theory of superconductivity was developed on the microscopic 
level 0. 

In the Standard Model, a number of striking results has been obtained using perturbation 
theory, among them the anomalous magnetic moment of the electron, asymptotic freedom 
in QCD and many others. In spite of these great successes, perturbation theory has some 
intrinsic problems: 

Perturbation theory is an expansion in powers of the gauge coupling. This is by def- 
inition only applicable to cases of small coupling. It is completely impossible to address 
strongly coupled problems in perturbation theory. Moreover, a perturbation expansion is 
an expansion about the trivial vacuum of the theory. But there exist objects in Yang-Mills 
theory, e. g. instantons, that are topologically distinct from the trivial vacuum. Such ob- 
jects can by no means be included in an expansion about the trivial vacuum. This can also 
be seen from the partition function: The instanton enters the partition function with the 
measure e~ s = e~ 8n I 9 , which has an essential singularity at zero coupling such that the 
Taylor series of this function about g = vanishes identically. Therefore, excitations of 
nontrivial topology are completely ignored in perturbation theory. 

For SU(N) Yang-Mills theory at finite temperature, further problems are known. It was 
shown that a perturbative calculation of the thermodynamical pressure can not be driven 
beyond order g 5 [I]. This is essentially due to the existence of weakly screened, soft magnetic 
modes which cause infrared instabilities. 

There is also a number of experimental and observational results in particle physics and 
cosmology the explanation of which is still open or being disputed. Examples are the non- 
detection of the Higgs particle at LEP, the existence of dark matter and dark energy in the 
universe, and a quadrupole signal in one of the power spectra of the cosmic microwave back- 
ground. A nonperturbative description of gauge theories may be helpful or even necessary 
in understanding these aspects. 

At present a large number of theories and approaches which try to include nonper- 
turbative aspects in field theories is being considered. An analytical and nonperturbative 
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approach to SU(2) and SU(3) Yang-Mills thermodynamics is presented in In the spirit 
of the Ginzburg-Landau theory of superconductivity, a macroscopic field is used to account 
for microscopic processes in an effective theory. More precisely, a composite adjoint Higgs 
field 0, which describes the BPS saturated and topologically nontrivial part of the ground 
state, is introduced. The Higgs field is generated at an asymptotically high tempera- 
ture by noninteracting calorons of topological charge one and trivial holonomy. The field 
is quantum mechanically and thermodynamically stabilized and can thus be used as a 
background for the topologically trivial sector of the theory. Interactions between trivial- 
holonomy calorons are included via a macroscopic pure-gauge configuration ajf which is a 
solution to the equation of motion for the topologically trivial sector in the presence of the 
background 0. 

As the modulus of the Higgs field decreases with temperature as \(f)\ ~ v/A 3 /T, where A is 
the Yang-Mills scale, the effects of topological defects die off at large temperature in a power- 
like fashion. Asymptotic freedom and infrared-ultraviolet decoupling of the fundamental 
theory, which are results obtained in perturbation theory at T = 0, are preserved. 

Thermo dynamical quantities are in this framework calculated as loop expansions about 
the nontrivial ground state consisting of the Higgs field and the pure-gauge configuration 
a b y . On tree level the excitations in the high temperature (or electric) phase are either 
massive thermal quasiparticles or massless 'photons'. The interactions between these quasi- 
particles appear to be very weak. The effective theory is both infrared- and ultraviolet-finite. 
The former is due to the existence of caloron- induced gauge boson masses (IR cutoff), the 
latter to constraints on loop momenta arising from the existence of the compositeness scale 
|0| (UV cutoff). 

The purpose of this thesis is to compute the dynamical generation of the macroscopic, 
composite field and to calculate the pressure of the Yang-Mills gas on two-loop level for 
the SU(2) case. The thesis is organized as follows: Chapter 2 reviews basic properties of 
Lie groups and pure Yang-Mills theory. The physics of some topological objects in gauge 
theories, namely the Abrikosov-Nielsen-Olesen vortex, the 't Hooft-Polyakov monopole, and 
instantons is sketched. The focus is on the latter. Chapter 3 first presents a brief outline 
of the physics of the electric phase according to the approach in [5j . The definition of the 
phase of the composite Higgs field in terms of a spatial and scale parameter average over 
an adjointly transforming two-point function is given and discussed. The average has to 
be evaluated on trivial holonomy caloron and anticaloron configurations. We discuss the 
uniqueness of the given definition and perform the evaluation. We show how under the 
assumption of an externally given scale the modulus of the field can be determined. The 
potential is deduced from the BPS equation. Chapter 4 contains the calculation of the 
two-loop contributions to the thermodynamical pressure of the SU(2) Yang Mills gas. We 
determine the contributing diagrams and state the Feynman rules. The computation of the 
diagrams is performed in the real time formalism of finite temperature field theory. The 
resulting integrals have to be evaluated numerically. We compare the two-loop contribution 
to the one-loop contribution and give an interpretation. Chapter 5 gives a short summary 
and an outlook on further research. The appendices contain technical details concerning 
the calculations in Chapters 3 and 4. 
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Chapter 2 



Basics of SU(N) Yang-Mills Theory 
and Solitonic Configurations 

2.1 Lie groups and Lie algebras 

For a continuous group G, the group elements may be parameterized by a set of continuously 
varying real parameters, a a (a = 1, 2, . . . , K). The total number of parameters K is called 
the order of the group. There are groups with compact parameter space such as SU(N), 
where K = N 2 — 1, and SO(N), where K = N ^~ 1 - > . There are also groups with non-compact 
parameter space, like the Poincare group in four dimensions, where K = 10. 
Lie group. A Lie group is a continuous group G where the set of parameters represents a 
different iable manifold. The latter is referred to as group manifold. The multiplication map 

GxG^G: (91,92)^91-92 (2.1) 
and the inverse map 

G - G : g - g' 1 (2.2) 

are differentiable. In terms of the parameters a and a' this means 

g(a) ■ g(a') = g(a") , (2.3) 

where a" is an analytic function of a and a' and similarly for the inverse map. As the 
dependence of the group elements on the parameters a a is analytic, any infinitesimal element 
g G G can be power expanded about the unit element of the group: 

g(a) = 1 + ia a T a + 0(a 2 ) . (2.4) 

The objects 

are called infinitesimal generators of the group G. For convenience, the unit element e has 
parameters a a = 0, g(0) = e. 
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Lie algebra. A vector space V together with a bilinear operation, the Lie-Bracket, given as 



V x V -> V : (X, Y) h-> [X, F] (2.6) 
and satisfying 

[X,X] = (2.7) 
and the Jacobi identity 

[X, [Y, Z]] + [Y, [Z, X}] + [Z, [X, Y]\ = (2.8) 
is called a Lie algebra. The property Eq. ()2.7|) implies the antisymmetry of the Lie bracket, 

[X,Y] = -[Y,X]. (2.9) 

The generators of a Lie group G always form a Lie algebra g. The dimension of the vector 
space V is equal to the number of generators (that is the order) of G. If G is a matrix group, 
then the Lie bracket is the usual matrix commutator and the Jacobi identity Eq. ()2.8|) is 
trivially fulfilled. 

Since the generators T a of the group provide a basis of the Lie algebra, the Lie bracket 
of two generators must again be a linear combination of generators, 

[T a ,T b }=if abc T c . (2.10) 

The symbols f abc are called structure constants. They can be chosen to be completely 
antisymmetric and real. Using the structure constants, the Jacobi identity for the generators 

[T a , [T b , T c ]] + [T b , [T c , T a ]} + [T c , [T a , T b }} = (2.11) 
can be phrased as 

jade jbcd _|_ jbde jcad _|_ jcde jabd g ^ \'2\ 

Many properties of a Lie group can be derived from the Lie algebra; e. g. the commutation 
relations Eq. ()2.10j) of a Lie algebra (which themselves often are called the Lie algebra) 
completely determine the multiplication law of the associated Lie group in the vicinity of 
the unit element. The connection between Lie algebra and Lie group is established through 
the exponential map, 

teg expteG, (2.13) 

which, in case of a matrix group, is the usual exponential map defined by the power series 

°° An 

ex p^ = E^- ( 2 - 14 ) 

n=0 

There are Lie groups that have the same Lie algebra but group manifolds of different global 
structure and topology. 

A Lie group G that contains no invariant sub-group (ideal) other than {e} and G itself 
is called simple. If G does not contain any Abelian invariant sub-group other than {e}, it 
is called semi-simple. The corresponding definitions apply to the Lie algebras. 
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For a Lie algebra, it is possible that some of the generators T a commute (i. e. their 
Lie-bracket is zero). These generators form the so-called Cartan sub-algebra, their number 
is known as the rank of the Lie algebra. For example SU(N) has rank N — 1; SO (2) and 
SO(3) both have rank 1. 

Representations . A linear representation R of a group G is a map which maps every group 
element g G G onto a linear transformation R(g) of a vector space V. This map has to 
respect the group multiplication, i. e. 

R(gi)R(g2) = R(gig 2 ) (2.15) 

and to map the unit element e onto the unit matrix, 

R(e) = l. (2.16) 

In other words, this map is a homomorphism. When a basis for V is chosen, R(g) are 
matrices. The dimension of the representation is by definition the dimension of the vector 
space V. In quantum mechanics, one is usually interested in finite dimensional unitary 
representations because unitarity is closely connected to probability density conservation. 

A representation is called reducible if it is possible to find a basis in which all represen- 
tation matrices have the form 

m=(^ (2.17) 

otherwise it is called irreducible. If additionally A(g) = for all g 6 G, the representation 
is called fully reducible. For a semi-simple group, all reducible representations are fully 
reducible (Weyl's theorem). 

Two representations R, R! of a group are said to be equivalent, if they only differ by a 
similarity transformation 

R>( g ) = S- 1 R(g)S VgeG (2.18) 



with a nonsingular matrix S. It can be shown that every representation of a compact group 
is equivalent to a unitary representation. 

To form a representation of the Lie algebra, the representation matrices have to respect 
the commutation relations Eq. (|2.1U|) of the Lie algebra. The elements of a matrix group 
can be viewed as linear transformations of lR n . Thus the group elements themselves form 
a linear representation of the group, the fundamental representation. The r generators of 
the Lie group G carry a representation of dimension r, the adjoint representation. The 
generator T a is mapped onto the mapping (denoted by the same symbol T a ) 

T a :g^Q; T h ^[T a ,T h \. (2.19) 



The representation matrices of the adjoint representation are given by the structure con- 
stants, 

R£i(T b ) =if abc . (2.20) 

As a direct consequence of the Jacobi identity Eq. (|2.8p . the adjoint representation ma- 
trices fulfill the Lie algebra commutation relations Eq. (|2.10J) . as requested. The adjoint 
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representation is always a real representation because the structure constants are real and 
antisymmetric. The adjoint representation of a simple Lie group is always irreducible. 
The group SU(N). SU(N) is the group of unitary N x N complex matrices with unit deter- 
minant, 

rfU=l, detl/ = 1. (2.21) 

It is a compact simple Lie group. Its Lie algebra has rank N — 1. The generators of SU(N) 
are the Hermitian and traceless N x N matrices 

H ] = H, tr# = 0. (2.22) 

The generators of SU(2) and SU(3) are usually taken to be the Pauli and Gell-Mann ma- 
trices, respectively. The adjoint representation of SU(N) has dimension N 2 — 1. 

2.2 SU(N) Yang-Mills theory 

We are considering a four- dimensional Minkowskian spacetime. An SU(N) Yang-Mills theory 
is governed by a Lagrangian which is invariant under any local SU(N) transformation. Every 
field in the theory has to transform under a unitary and finite dimensional representation of 
SU(N). A Yang-Mills theory containing only gauge fields but no matter fields is often called 
pure. 

The ordinary derivative cannot be used to construct gauge invariant quantities. So, 
in order to be able to include derivative terms (i. e. kinetic terms for the fields) in the 
Lagrangian, the covariant derivative and a gauge field have to be introduced 1 , 

= <9 M - ieA„ . (2.23) 

The gauge field A^ can be expanded in terms of the Hermitian generators of SU(N), 



A, = A a ^ , (2.24) 



where the SU(N) generators A a are normalized such that trA a A b = 25 a b- When applying 
the covariant derivative to a matter field, the gauge field A^ is understood to act in the 
representation of the matter field. That is for a fundamental field ip 



X a 

and for an adjoint field 



D^^d^-ieA^—cp (2.25) 



(D^) a = d^ a + ef abc A b ^ c . (2.26) 

The covariant derivative is constructed such that the covariant derivative of a field has 
exactly the same transformation law as the field itself. To satisfy this request, the gauge 
field has to transform under an SU(N) gauge transformation U(x) according to 

A^x) -> U(x) A^x) U\x) + ~ U(x) d» U\x) . (2.27) 



1 One can conventionally absorb the gauge coupling e in the definition of the gauge field and write 
Dft = — iA^ and F^ v = d^A v — dyA^ — ifA^AJ etc. This notation is convenient for considering 
nonperturbative aspects. We will use it for working with instantons in Sees. 12.51 12.61 and in Chapter 3. 
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Viewing the gauge group (being a Lie group) as an analytic manifold, the gauge field 
is a connection on this manifold. The curvature of the manifold (in differential geometry) 
corresponds to the field strength tensor (in field theory), namely 

Ffiv = % ~ [D„, D v \ = d„A v - d u A^ - xe [A„ A u ] . (2.28) 

In contrast to Abelian theories, the field strength in nonabelian theories is not a gauge 
invariant quantity but transforms under the adjoint representation as 

F^{x) -> U(x) F^(x) U\x) . (2.29) 

The field strength can also be written in matrix notation, 

Y 

with the components 



F,u = F; v - , (2.30) 



F% = d,>AZ - d v Al + ef ahc A\Al , (2.31) 

where f abc are the structure constants of the gauge group. The kinetic term for the gauge 
field A^ in the Lagrangian is 

C = -itrF^F"" = —F^F^. (2.32) 

If one demands Lorentz invariance, gauge symmetry, renormalizability and CP-invariance, 
no further terms are allowed in pure Yang-Mills theory. Relaxing the demand for CP- 
invariance, an additional term proportional to F^ v F^ ua may be added to the Lagrangian in 
Eq. ()2.32|) . Here F^ u = \e nuapF a!3 denotes the dual field strength. 

In particular a mass term for the gauge field is forbidden because it has the gauge variant 
form m 2 tr A^A^. Nevertheless, gauge bosons can acquire mass by dynamical symmetry 
breaking which is manifested by the Higgs mechanism. 

The commutator term in Eq. ()2.28|) vanishes for Abelian gauge groups, so the Lagrangian 
in this case is quadratic in the gauge field, the equations of motion are linear in the gauge 
field, and hence there is no self interaction. If, in contrast, the gauge group is nonabelian 
(as is SU(N)), this is no longer true. The Lagrangian contains terms cubic and quartic in 
the gauge field (besides quadratic terms) and thus allows for three- and four-gauge boson 
vertices. Because of this a pure SU(N) Yang-Mills theory is interacting. 

The equations of motion for the field A^ are derived via the minimal action principle. 
For SU(N) pure Yang-Mills theory, they are 

D^ u = (2.33) 

or in components 

d^F^ + ef abc A b ^F^ = . (2.34) 

The right hand side of the equation of motion is zero because of the absence of external 
sources, i. e. charged matter. Nevertheless, this is an interacting theory since the gauge field 
couples to itself due to the nonlinearity in the field tensor. 

On the classical level the Lagrangian Eq. ()2.32j) does not contain any dimensionful pa- 
rameter: The gauge coupling e is dimensionless, and gauge boson masses are forbidden. 
This is the reason for the invariance of the action under a rescaling of fields and spacetime 
arguments. On the quantum level, however, a mass scale comes into existence due to the 
mechanism of dimensional transmutation [Zj. 
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2.3 The Abrikosov-Nielsen-Olesen vortex line 

Consider a theory with a U(l) gauge field and a scalar Higgs field <fi defined by the 
Lagrangian 

C = -~F^F» V + (D,<p)(Dy) - m 2 0*0 - A(0*0) 2 . (2.35) 

If m 2 < 0, the gauge symmetry is spontaneously broken, and the Higgs field acquires a 
vacuum expectation value 



-m? 
~2\~ 



(2.36) 



The Higgs Lagrangian Eq. ()2.35|) exhibits static string-like solitonic solutions, so-called vor- 
tices. They were discovered by Nielsen and Olesen in 1973 [8 . Consider a field configuration 
with cylindrical symmetry and the asymptotic behavior 



|vac e me (r^oo) (2.37) 
for the Higgs field, and 

= -^{nO) (r -> oo) (2.38) 

for the gauge field. In cylindrical coordinates z, r = \x\ and 9, the latter reads 



A z = 
A r = 

A = (r -> oo) 

er 



(2.39) 



Because of the demand for single valued fields, n has to be an integer. In the asymptotic, 
An is a pure gauge and thus the field strength vanishes, 

F^ u = (r -»■ oo) . (2.40) 

The Higgs field <fi is covariantly constant, 

D r <\) = and D e <p = (r -»• oo) , (2.41) 

and the potential V = m 2 <j)*<f) + A(0*0) 2 evaluates to zero. So the energy density at r — > oo 
is TC = —C = 0, and static solutions with finite energy and the above boundary conditions 
can exist in principal. The equations of motion obtained from the Langrangian Eq. ()2.35|) 
are 

D^(D^) = -m 2 - 2A0|0| 2 

^e(0^0*-0*^0)+2e%|0| 2 = ^F^. l ' J 

One can check that the equations of motion allow for configurations with the asymptotic 
behavior Eqs. ()2.37j) and (|2.38jl . To find a solution to the equations of motion which satisfies 
the above boundary conditions, one makes the ansatz 

A z (r) =0 

A r (r) = (2.43) 
A B (r) = A(r) 



for the gauge field, and 

<P = X(r)e md (2.44) 

with 

X(r) — -> and x(r) ► Mvac (2.45) 

for the Higgs field. The magnetic field B will have only a z-component, 

B t = ~[rA(r)]. (2.46) 
r dr 

Inserting the above ansatz into the equations of motion yields a system of differential equa- 
tions for the functions A(r) and x{ r )i 



1 d ( dx 

r 



r dr \ dr 
d_ 

dr \rdr 



n \ ^ 

eA) + m 2 + 2Ax 2 

r / 



(~(rA))-2e (- + eA) X 2 = . 
\r dr J \e / 



x = o 

(2.47) 



No exact solutions to these equations are known. The asymptotic behavior of the gauge 
field and the magnetic field has been deduced by Nielsen and Olesen as 

A = —— - -^(|e||0Ur) ► - - ( * V /2 e-NI0U, + . . . (2 . 48) 

er e ' r^<x> er e \2\e\ |0| vac rj 



and 

2 

r^oo e \ 2 e r 



B z = cxK o (|e||0| vac r) ► £ ( Bf!p ) C -WW— ■ + . . . , (2.49) 



where Kq and denote modified Bessel functions, and c is a constant of integration. 

The line integral over around a circle S\ at infinity yields the magnetic flux through 
the surface enclosed, 

$ = J B da = j Apjdx^ = j A e rd9 = -~n , (2.50) 

The magnetic flux is quantized: it appears only in multiples of the flux quantum — . 

The vortex line owes its existence to the topological structure of the gauge manifold. 
The boundary condition Eq. (|2.37|) defines a mapping of the boundary S\ in physical space 
onto the group manifold of U(l), which again is S\. There are infinitely many classes of 
such mappings which can not be continuously deformed into one another, 

7r 1 (Z7(l))=7r 1 (5' 1 ) = Z. (2.51) 

A gauge theory, where the gauge group G has iti(G) = 0, does not exhibit vortex lines. This 
is the case for SU(2), for example. 

The Lagrangian Eq. ()2.35j) is the relativistic generalization of the Landau-Ginzburg free 
energy of a type II superconductor. The electromagnetic field A^ interacts with the bosonic 
field 0, which describes the Cooper pairs. In the superconducting phase the photon becomes 
massive, and if an external magnetic field can enter the superconducting material (i. e. a 
type II superconductor), it does so only in the form of (quantized) Abrikosov flux tubes. 
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2.4 The 't Hooft-Polyakov monopole 

Magnetic monopoles in gauge theories were first considered in 1974 by 't Hooft jU] and 
Polyakov JU] i n the context of an SU(2) Yang-Mills theory with an isovector Higgs field 4> a . 
The Langrangian of such a theory is 

I 1 TP 2 

with the parameter m 2 chosen negative such that the Higgs field has a non-zero vacuum 
expectation value F with 

F 2 = --. (2.53) 

Now, consider static solutions with the asymptotic behavior 

A a i =-e iah ^ i (r^oo) 

er 2 (2.54) 

A a = 

for the SU(2) gauge field and 

<f) a = F— (r -> oo) (2.55) 
r 

for the Higgs field. Here r is the norm of the spatial vector, r = \x |. Note that in Eqs. (|2~54|l 
and 1)2.55)1 space and isospace indices are mixed. Fields of the form Eq. ()2.55j) are known as 
" hedgehogs" . 

By the presence of a nonvanishing vacuum expectation value of the Higgs triplet the 
SU(2) gauge symmetry is broken. The field strength corresponding to the unbroken U(l) 
subgroup is the 't Hooft tensor 

F, u = ^> a F% - -^e^D^D^) . (2.56) 

Upon inserting 



A 1 ^ 2 — A 3 ^ — A^ 

6 1 ' 2 = 6 3 = F , 



(2.57) 



it reduces to the usual definition of the electro-magnetic field tensor. Inserting the nontrivial 
asymptotic conditions Eqs. ()2.54)) and ()2.55)) the 't Hooft tensor evaluates to 

F^ = -—.e^x* (r -> oo) . (2.58) 
This corresponds to the radial magnetic field of magnetic point charge, 

B a = ^ (r -> oo) (2.59) 
with a total flux or magnetic charge 

/ 47T 

q mag = $ = <p da a B a = — . (2.60) 
Js 2 e 
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In [0] it is shown that there exist configurations with the requested asymptotic behavior that 
are smooth and hence have finite energy for all A and m 2 < 0. In fHJEi] they are explicitly 
given for the so-called BPS limit A — > 0. In the BPS limit, the mass of the monopole is 

M m = ^M w , (2.61) 

where M\y = eF = is the vector boson mass. For general A, the mass of the monopole 
is larger but still of the same order. 

The magnetic charge can also be expressed as 

Qm ag = ^ f d 3 xK° (2.62) 

with the current 

K" = — d wpa e abc d u 4> a dj b dj c 

2e p (2.63) 

= d F 

where (j) a = ^ . This current is identically conserved, d^K^ = 0; it is not a Noether current 
corresponding to some symmetry of the Lagrangian, but of a topological nature. Inserting 
Eq. ()2.63|) into Eq. ()2.62|) and applying Gauss' theorem yields 

q m = -3— / d(Ti e ijk e abc 4> a dj4> b d k (j) c , (2.64) 



we 



S-2 



where the integral has to be performed over an S2 with infinite radius. From Eq. (|2.64j) it 
can be seen that the total magnetic flux is completely carried by the Higgs field. Moreover, 
Eq. 1)2.64)1 states that the magnetic charge does only depend on the asymptotic behavior 
of the fields. The existence of configurations with non-zero magnetic charge is due to the 
possibility to demand nontrivial boundary conditions for the fields, that means: There are 
maps from the surface of space (S2) onto the manifold S2 (corresponding to rotations in 
isospace) which can not be continuously deformed into the constant map. Mathematically 
speaking, the second homotopy group of S2 is nontrivial, 

7r 2 (S 2 )=Z. (2.65) 

Conservation of magnetic monopole charge is due to this topological argument, namely the 
transition between nonhomotopic gauge configurations needs infinite energy. 



2.5 Instant ons at zero temperature 

Introductory material on instantons and other topological objects in gauge theories is pre- 
sented in |TT|. Reviews on instanton physics are jT2| IT3]. 

SU(N) vacuum. The (Minkowskian) vacuum of an SU(N) Yang-Mills theory is infinitely 
degenerate. There is not only the trivial vacuum A^(x) = 0, but every field configuration 
of the form 

A^x) = iU(x)d ll U\x), (2.66) 
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Figure 2.1: Left-hand side: In the potential V(x) the transition from A to B is a tunneling 
process and hence classically impossible. Right-hand side: In the reversed potential —V(x) 
the transition from A to B is classically possible along the indicated trajectory because 
allowed and forbidden region have been interchanged. 



where U(x) is an SU(N) matrix, differs from = only by a gauge transformation and 
hence has zero field strength and energy density, as well. Configurations as in Eq. ()2.66|) are 
referred to as pure gauge. These pure gauge configurations fall into topologically distinct 
classes, and hence they cannot be smoothly connected. The configurations Eq. (|2. fifty are 
classified by a winding number, the Pontryagin index, defined as 



n w 



24tt 2 



dPxe* jk tT(rfd i U){rfd j U){Uid k U) . (2.67) 



The quantity nw is an invariant under smooth deformations. It can also be expressed in 
terms of the gauge field as 



n w 



167T 



J d 3 xe^ k UtdjAt + ^f abc A1A)A% \ . (2.68) 



It is not possible to find a solution to the equation of motion with finite energy that connects 
two vacua with different winding number. Configurations corresponding to a transition 
between different vacua of the theory via a tunneling process are called instantons. They 
have to be considered in the framework of Euclidean field theory. 

Tunneling solutions. Tunneling configurations are solutions to the Euclidean equations of 
motion (i.e. they minimize the action of the theory, formulated in a Euclidean spacetime). 
This can be motivated by an example from the quantum mechanics of a point particle. 

Imagine a point particle with kinetic energy T moving in a time independent potential 
V. A semiclassical approximation (that is an expansion around the classical path) will not 
include any tunneling process because there are no classical tunneling paths. Now letting 

E-V^V-E (2.69) 

clearly interchanges allowed and forbidden region, so that a classical path exists, see Fig. 12.11 
An expansion around this classical path yields the WKB approximation for the tunneling 
process. But E — V = T, and T is the square of a time derivative. Hence the sign change 
Eq. ()2.69|) is equivalent to the Euclidean continuation, 

t -> ir (r real) . (2.70) 
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So tunneling processes are solutions to the Euclidean equations of motion. This is true in 
quantum field theory as well. 

Finite action and boundary conditions. In order to get finite action configurations, one 
has to demand that the field strength F^ u approaches zero at the boundary of spacetime 
faster than 1/x 2 . This can not only be fulfilled trivially (by approaching = 0), but 
also nontrivially by demanding that the gauge field is pure gauge on the boundary of 
Euclidean spacetime (which is S3), 

A tl = iU(x)d fl U^(x) (\x\->oo). (2.71) 

This condition defines a mapping from the surface S3 to the group space of SU(N). The 
equivalence classes of homotopic mappings from S3 into a manifold X, f : S3 — > X, form 
the so-called third homotopy group of the manifold, 7r 3 (X). So the finite action field config- 
urations fall into topologically distinguished classes if and only if the third homotopy group 
of the gauge group manifold is nontrivial. The homotopy group n 3 of most gauge groups is 
known, see for example [14]. The case of basic interest is SU(2), where 

7r 3 (SU(2)) = 7r 3 (S 3 ) = Z. (2.72) 

The integer attached to a given mapping / : S3 — > S3 indicates how often the gauge group 
manifold S3 is "wrapped" around the spacetime-S3 and hence is referred to as winding 
number. It is also known as the Brouwer degree of the mapping /. 

Pontryagin index and topological charge. For a given field configuration the four dimen- 
sional Pontryagin index is defined as 

Q = -zkzl ( 2 - 73 ) 



32vr 2 J " u <* u 

The integrand of Eq. (|2.73jl . the Pontryagin density, can also be written as a total divergence, 

^ Fa J% = d ^ > ( 2 - 74 ) 
where is the Chern-Simons current, 

= 16^ £ ^ ( A " d e A " + \f abc KA h p A^ . (2.75) 

In contrast to the Pontryagin density, the Chern-Simons current is a gauge variant quantity. 
If the integrand is nonsingular, the volume integral in Eq. ()2.73j) can be converted to a 
surface integral by means of Gauss' theorem, 

Q = J d 4 xd^K„ = J da^K^. (2.76) 

Although this surface integral has to be evaluated on a sphere S3 with infinite radius, it is 
not necessarily zero: If the gauge field A^ falls off less than l/|x|, then the Pontryagin index 
will be non-zero. 

From Eq. (|2.76J) it is obvious that (for regular gauge fields) the Pontryagin index is 
solely determined by the asymptotic behavior of the gauge field or in other words by the 
boundary conditions. Thus the existence of fields with non-zero Pontryagin index is due 
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to the possibility of demanding nontrivial boundary conditions, as already indicated in the 
context of Eq. (|2.72j) . The Pontryagin index is a topological charge because its conservation 
does not follow from a continuous symmetry of the Lagrangian but is solely due to the 
stability of the boundary conditions under continuous perturbations. 

Bogomolnyi bound and self-duality. For the construction of a configuration of minimal 
Euclidean action connecting two vacua with different winding (Chern-Simons) number, the 
Yang-Mills action is written in form of a Bogomolnyi decomposition, 

s = At / J/: 1 .. 



Ag 2 J ^ ^ 

i r / i \ ( 2 - 77 ) 



where 



F^v — -E^apFap (2.78) 

is the dual field strength tensor. The second term in Eq. ()2.77|) is a square and hence always 
positive; the first term is the Pontryagin index Q defined as in Eq. ()2.73|) . a topological 
invariant. Hence any self-dual (or anti self-dual) field configuration, i.e. 

Fpv = ±F IIV , (2.79) 
is a minimum of the Euclidean action (in a given topological sector). It has the action 

S=%\Q\. (2.80) 

r 

In the dual field strength the roles of electric and magnetic fields are interchanged as 
opposed to F^y. So, (anti) self-duality can (at least in temporal gauge) also be characterized 
by 

E a = ±B a , (2.81) 

where E a and B a are the color electric and magnetic fields respectively. From the decompo- 
sition Eq. ()2.77|) it is clear that a self-dual gauge field is a (local) minimum of the Yang-Mills 
action, and indeed self-duality and the Bianchi identity imply that the equations of motion 
are satisfied, 

D fl F llu = ±D fl F llu = 0. (2.82) 

In contrast to the equation of motion ()2.82|) . which is a second-order differential equation, 
the self-duality equation (j2.79j) is first order. The Bogomolnyi decomposition Eq. ()2.77j) 
gives a lower bound for the action in a given topological sector. Configurations saturating 
this bound (i.e. self-dual configurations) are also referred to as BPS-saturated, [To"! ITd']. 
It can be shown that the energy-momentum tensor vanishes identically on BPS-saturated 
fields. 

BPST instanton. Here we consider the gauge group SU(2). In the boundary condition 

-> iUdjfl (2.83) 
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which is necessary to achieve a configuration of finite Euclidean action (see above), one 
chooses 

X4+^A flsl^oo) (2.84) 
\x\ 

where |x| = y/x^x^, and A, the are Pauli matrices. Therefore the gauge potential has to 
have the asymptotic form 

A l^H»% (M->°°). ( 2 - 85 ) 

where the 't Hooft symbols 77" and 77" are defined as 

V/iu = e a^u + 8<x)jfiu4 — 8 a v5fi4 (2 86) 

This leads to the ansatz 

^ = 2 ^^' ( 2 - 87 ) 

where the scalar function / has to satisfy the boundary condition / — > 1 for \x\ — > oo. 
Inserting this ansatz into the self-duality equation (|2.79jl yields the differential equation 

/(l-/)-xY = 0, (2.88) 

which is solved by 

where p is a constant of integration. This gauge field configuration is the Belavin-Polyakov- 
Schwartz-Tyupkin instanton [TJj 

A » = 2 ^^f? ' (2 ' 90) 

For this instanton solution, the Chern-Simons current Eq. (|2.75J) has the asymptotic behavior 



K »-+ww (N ^ oo) ' (2 - 91) 

and inserting this into Eq. ()2.76|) yields Q = 1 for the topological charge of the BPST 
instanton. Replacing 77" with 77° yields a solution with Q = —1. The corresponding field 
strength is 

P 2 

F W = ( x 2 _)_ ^2)2 ' (2.92) 

Although the instanton gauge potential only falls of as 1/x, the field strength Fj} v falls of 
as 1/x 4 . So instantons are well localized in (Euclidean) time and in space. 
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The solution Eq. (J2.90)) can be generalized by shifting the instanton center to an arbitrary 
point Zp, 



and changing the color orientation of the instanton by a global gauge rotation (containing 
three parameters). Furthermore the BPST instanton contains an arbitrary parameter p, 
which gives the size of the instanton. These make up eight parameters in the generalized 
BPST instanton. Since the action of the configuration Eq. (j2.9Hj) . namely 

*=£|«| = £, (2.94) 

is independent of these eight parameters, they correspond to zero modes. The gauge in 
which the BPST instanton has the form (J2.90)) is referred to as regular gauge. 
Singular gauge. The BPST instanton Eq. ()2.90|) has a singularity at infinity. This singularity 
can by means of the singular gauge transformation 

\x\ 

be shifted from infinity to the origin, such that the instanton in so-called singular gauge has 
the gauge potential 



A » = = -%» ln (l + ) • (2-96) 



The singularity of the gauge field at x = is not physical; field strength, action density 
and topological charge density are smooth. The field strength of an instanton in singular 
gauge is 



w = -prW («• " ^ir " 2C - £ f i ) • (2 - 97) 



In regular gauge the integral over the topological charge density picked up only contributions 
at |x| — > oo, cf. Eq. (|2.76j) . In contrast to that, in singular gauge not only the action density, 
but also the topological charge density is localized at the center of the instanton. It is this 
property that allows for the construction of multi instantons in singular gauge. 
Multi instanton solutions. The 't Hooft ansatz for the gauge potential 

A»(x) = -^^lnn(x) (2.98) 

is a generalization of the BPST instanton in singular gauge, Eq. ()2.96|) . The demand for 
self-duality translates into a condition for the scalar pre-potential H(x), it has to fulfill the 
Laplace equation 

n- 1 (x)-^n(x) = o. (2.99) 

Hence there is (only in singular gauge) some kind of superposition principle: Superposing the 
pre-potentials of single instantons with different centers and scale parameters again yields a 
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self-dual solution to the equations of motion, a multi instanton. Therefore the pre-potential 
n(x) has the general form 



K 




U(x) = i + J2 



(x - z n ) 2 ' 



(2.100) 



n=l 



The configuration Eq. ()2.98j) with the pre-potential Eq. ()2.100|) inserted has Pontryagin index 
Q = K. Replacing rf^ v in the 't Hooft ansatz with 77" yields a solution with Q = —K. The 
't Hooft ansatz in particular allows for the construction of periodic instantons. 
Collective coordinates. The collective coordinates of the SU(2) instanton with topological 
charge \Q\ = 1 are the instanton size p, the instanton position z and three parameters 
determining the color orientation, making up 8 parameters in total. The action on an 
instanton is independent of all theses parameters; i. e. they are moduli of the instanton 
solution and correspond to zero modes. 

A multi instanton of topological charge K is the superposition of K single instantons and 
accordingly has 8K parameters, [T^j. Note that the multi instanton solution as displayed in 
Eqs. (|2.98J) and (|2. 100J) does not cover all of this moduli space: the K instantons are chosen 
to have the same color orientation, which must not necessarily be the case. The general 
SU(2) multi instanton solution has been constructed by Atiyah et. al. [TH] . 
Gauge groups other than SU(2). For the case of a simple gauge group G, basically no 
instantons other than in SU(2) arise, but the various possible embeddings of SU(2) in G 
have to be taken into account. The general SU(N) multi- instanton with Pontryagin index 
Q has 4NQ parameters, pjlTS]. 

Tunneling interpretation. Assuming that the (regular) gauge potential falls of rapidly at 
spatial infinity, one can write the topological charge Q from Eq. (|2.7fij) as 



= n w (r = 00) - %(r = -00) , 

where the fact that the zero component of the Chern-Simons current Eq. ()2.75j) indeed is 
equal to the integrand of Eq. ()2.68|) was used. That means an instanton with Pontryagin 
index Q 7^ connects topologically different vacua with the according difference in winding 
number via tunneling. 

2.6 Instantons at finite temperature 

Finite temperature. The equilibrium thermodynamics of a quantum field theory at finite 
temperature usually is considered in the framework of Euclidean field theory. The finite 
temperature T corresponds to imaginary time being compactified on a circle with circum- 
ference /3 — h. The physical fields have to be periodic in Euclidean time. For a gauge 
field the periodicity condition is less strict; only periodicity up to a gauge transformation is 



Q 




(2.101) 
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demanded, 



Af,((3, x) = fi(0, x) A^O, x) Q\0, x) + -{1(0, x) Q\0, x) (2.102) 

with Q e SU(N). We now have to discern time and space and write = (r, x). 
Periodic instanton. To find an instanton solution which is periodic in Euclidean time, one 
starts from the 't Hooft ansatz Eq. ()2.98|) with the pre-potential Eq. ()2.100|) describing a 
multi instanton built out of K instantons with scales p n and centered at z n respectively. 
All of them have the same color orientation. Taking the sum over infinitely many instanton 
pre-potentials centered at (n/3, 0), n e Z with the same scale parameter p yields a periodic 
instanton (or caloron) centered at (0,0), 

oo 2 

n(r,*) = i+ E (fz^iy - < 2103 > 

n=— oo v ii/ 

The sum has been performed by Harrington and Shepard, see [TH] . 



7rn2 sinh 

n ^ 1 + i cosh ¥r -eosy, ' (2 ' 104) 



where r = \x\. The caloron of Eq. ()2.104j) still has the action S = the action is, 
on the classical level, independent of temperature. By superposing (the pre-potentials of) 
single-calorons a caloron solution of higher topological charge can be constructed. As in 
the case of the multi-instanton, this is only possible in singular gauge. Calorons of higher 
topological charge also possess more moduli, some of which (namely the core sizes and the 
core distances) are dimensionful. 

Polyakov loop. The Polyakov loop is a time-like Wilson loop 



P(a;) = Pexp 



l/T 

ig I &t A±{t,x) 



o 



(2.105) 



where V denotes the path ordering operation. It has to be evaluated in periodic gauge, 
Aplr + P, x) = Ah(t, x). The Polyakov loop is only defined at finite temperature. At spatial 
infinity, the Polyakov loop does not depend on x/r any longer. Its value P(|cc| — > oo) is a 
topological invariant. 

Calorons are classified according to the eigenvalues of their Polyakov loop at spatial 
infinity. By definition, trivial holonomy means for SU(2) calorons that P(|aj| — ► oo) = 
±1. Calorons with non-unity eigenvalue of P at spatial infinity are said to have nontrivial 
holonomy, accordingly. The caloron solution from Eq. ()2.104j) has Polyakov loop P = 1. 
More general solutions have been constructed by Nahm [20], Lee, Lu j^I], and Kraan, van 
Baal [221 

Monopole constituents. In 1998 it has been shown independently by Lee and Lu |2I] and 
Kraan and van Baal [221 ESj that SU(N) calorons of nontrivial holonomy contain constituent 
BPS magnetic monopoles. These monopoles are subject to an attractive interaction in the 
case of small holonomy (i. e. close to trivial holonomy), or to a repulsive potential for large 
holonomy (i.e. far from trivial holonomy), see [21]. The latter can lead to dissociation of 
the caloron into a monopole-antimonopole pair. 
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Chapter 3 



Composite Adjoint Higgs Field in 
SU(2) Yang-Mills Theory 

In [5] an analytical and nonperturbative approach to the thermodynamics of SU(2) and 
SU(3) Yang-Mills theory in four dimensions is developed. An essential ingredient is the 
existence of an adjoint Higgs field composed of trivial holonomy calorons. 

Sec. 13 .11 gives a brief outline of this approach as presented in p]. We confine ourselves to 
the case of SU(2) and its so-called electric phase. We will give and evaluate a microscopic 
definition of the composite field m Sec. 13 .2\ and briefly comment on the nonexistence of 
alternative possibilities for its definition. 



3.1 An outline of the approach in [3| 

Decomposition of gauge fields. Each SU(2) gauge field appearing in the partition function 
of the fundamental theory is uniquely decomposed into a topologically nontrivial and BPS 
saturated part A l ° p represented by calorons, and a topologically trivial remainder a M , 

A, = A^ + a,. (3.1) 

At a large temperature, calorons with topological charge one and trivial-holonomy are as- 
sumed to generate a macroscopic adjoint scalar field 0. Interactions between the trivial 
holonomy calorons via topologically trivial fluctuations are not included in the field but 
are accounted for at a later stage by means of a pure-gauge background ajf. The field has 
to have the following properties: 

1. It describes (part of) the ground state of a thermal system, so its (gauge invariant) 
modulus |0| must be independent of space and time. 

2. The action of a classical caloron is independent of temperature. So no explicit T 
dependence may arise in 0's definition. 

3. Calorons are BPS saturated (or self-dual) solutions to the Yang-Mills equations of 
motion in four- dimensional Euclidean spacetime (with time r compactified on a circle). 
In particular, their energy-momentum tensor is precisely zero. The adjoint scalar field 0, 
being composed of noninteracting calorons, inherits this property, which in turn is expressed 
through a BPS equation for 0's time dependence, 

«9 T = V {1/2) . (3.2) 
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Here, is a 'square-root' of the potential V(4>) = tiiV^^^yV^ 1 ^ 2 ^, which governs the 

dynamics of 0. The potential V(<f>) is determined by the demand for BPS saturation. 

In Sec. 13.21 we will give a microscopic definition for and after evaluation see that this 
field is (up to a global gauge rotation) of the form 



A 3 

(r) = V 2^T Al ex P(- 2 ™ TA 3^) • (3-3) 

The potential is found to be 

V(4>) = A 6 tr cj)~ 2 = AttTA 3 . (3.4) 

Here, A is a fixed mass scale generated by dimensional transmutation. The modulus \(f)\ falls 
off as ;^=. This dependence shows that the nontrivial-topology sector is strongly suppressed 
at large temperature. 

Effective theory. Minimally coupling of <fi to the (up to now not included) topologically 
trivial fluctuations a M results in the effective action for the electric phase, 

S E = ^ dr j d 3 x QtrG^G^ + trV^V^ + A 6 0~ 2 ) , (3.5) 

where 

G% = ~ dual - ee abc alal (3.6) 

is the field strength of topologically trivial fluctuations a M , 

£> M = dp + ie [(f), a M ] (3.7) 

is the covariant derivative, and e denotes the effective gauge coupling. There is no reason 
why the effective gauge coupling constant should be equal to the coupling constant of the 
fundamental theory. 

The field <fi is seen to be quantum mechanically and statistically inert: The mass associ- 
ated with its excitations is much larger than both temperature and the scale of admissible 
quantum fluctuations of \(j)\ itself, 

$£ = IV and ^ = 3A*. (3.8) 

where A = is the dimensionless temperature. In the electric phase of the theory, A > 
11.65 (see below), so that 3A 3 indeed is a large number. The quantum mechanically and 
thermodynamically stabilized field can now be taken as a background for the equation of 
motion governing topologically trivial fluctuations, 

V fl G flv = 2ie[<t>,V v <f>). (3.9) 

This solution is to be part of the ground-state description and hence must not break the 
rotational invariance of the system. So it has to be a pure gauge. A pure-gauge solution to 
Eq. (pSl reads 



a b 9 = 5^T\ 3 . (3.10) 
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It takes into account holonomy changing interactions between calorons mediated by topo- 
logically trivial fluctuations. Moreover, we have 



2V> = 



(3.11) 



on the ground-state configuration 0, eft. As a consequence, the action density in Eq. ()3.5|) 
evaluates to V(<f>) on the ground-state. This corresponds to a ground-state energy density 
Pg.s. = V(<j>) = 47rA 3 T and a ground-state pressure P g , s , = —V(4>) = —4tiA 3 T. 

The inert scalar and the pure-gauge solution erf form the ground state about which 
loop expansions are performed. Therefore, the topologically trivial sector a M is split into the 
pure-gauge ground-state part ajf and fluctuations 5a M about this background. The adjoint 
scalar renders some of the fluctuations 5a M massive through the adjoint Higgs mechanism. 
In the case of SU(2), the gauge symmetry is dynamically broken 1 to U(l). Two of the gauge 
bosons acquire the temperature dependent mass 



while the third gauge boson is left massless 2 . The former are referred to as tree- level 
heavy (TLH) modes, the latter as tree-level massless (TLM) modes. Since the masses of 
TLH excitations are temperature dependent, they are thermal quasiparticle fluctuations. 
The masses induced by the Higgs mechanism provide infrared cutoffs, so that no infrared 
divergences occur in loop expansions of thermodynamical quantities. 

Compositeness constraints. The existence of the compositeness scale |0| present in the 
vacuum additionally implies cutoffs for the momenta propagating as vacuum fluctuations: 
1. The offshellness of a quantum fluctuation must not exceed the scale |0|, 



where p is the four-momentum of the quantum fluctuation and m its mass. 

2. The center-of-mass energy flowing into or out of a four- vertex must not exceed the 
scale |0|, 



where pi are the ingoing momenta. 

The violation of any of these two conditions would immediately imply that the fluctu- 
ations generated out of the vacuum are able to destroy this vacuum, which is impossible. 

1 With the introduction of a composite field which breaks the gauge symmetry partially this approach 
conceptually resembles the Landau-Ginzburg-Abrikosov theory of superconductivity Q][2]- 

2 The calculation of the mass spectrum as well as other explicit calculations (such as the thermodynamical 
pressure in Chapter [3J are carried out after a transformation to unitary gauge, where a^f = and hence 




(3.12) 



\p 2 -m 2 \ < |0| 2 



(3.13) 




(3.14) 



The scalar field in this gauge takes the form 




It can be shown that the employed gauge transformation is admissible, see 
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Both conditions Eqs. ()3.13j) and (J3.14)) have been formulated for a Minkowskian signature 
here, but can be rotated to the Euclidean. 

By virtue of the cutoffs as they emerge above, the phase space for vacuum diagrams 
is strongly restricted. As a consequence of this, the interactions between quasiparticle 
fluctuations are very weak and (for sufficiently large temperatures) generate only tiny higher 
loop corrections to thermodynamical quantities. Moreover, ultraviolet divergences do not 
occur and the usual renormalization procedure is not needed in the effective theory. 
Thermodynamical self-consistency. The thermodynamical quantities pressure, energy den- 
sity, entropy density etc. as obtained from the fundamental SU(N) Yang-Mills theory are 
related by Legendre transforms. These relations hold in general, so their validity in the effec- 
tive theory described by the action Eq. ()3.5j) has to be arranged for by imposing a condition 
of thermodynamical self-consistency. This condition has to assure that the T-derivatives of 
the TLH masses and the ground-state pressure cancel one another, so that only the explicit 
T-dependence arising from the Boltzmann weight enters T-derivatives of thermodynamical 
quantities. In particular, the relation 

dP 

P = T--P (3.15) 

translates into 
dP 

- = 0, (3.16) 

where 



, 3/9 / A 3 m , . 

a = 2ne\- 3/2 = e\ —- = — 3.17 
V 2ttT 3 2T v ; 

is a dimensionless measure for the quasiparticle mass. Eq. (jH.lfij) results into an evolution 
equation A (a) for temperature as a function of the tree-level gauge boson mass. The evo- 
lution has two fixed points, namely a = and a = oo. They signal the existence of both 
a lowest and a highest attainable temperature in the electric phase. These are denoted as 
X c = X(a = oo) and Ap = A (a = 0) respectively. The evolution is subject to the initial con- 
dition A(a = 0) = Ap. The low-temperature behavior of A(a) is practically independent of 
Xp as long as Ap is sufficiently large: The temperature at which the field emerges does not 
influence the low-temperature physics. This is a signal of ultraviolet- infrared decoupling. 

The evolution equation A(a) can be inverted to yield an evolution e(A) for the effective 
gauge coupling as a function of temperature. For SU(2) and SU(3) this is shown in Fig. 13. II 
The critical temperature is A c = 11.65 for SU(2), and the evolution exhibits a plateau, where 
the effective gauge coupling has the value e = 5.1. 

Phase transition to magnetic phase. At the critical temperature A c , the theory undergoes a 
second order phase transition from the electric to the so-called magnetic phase. The effective 
gauge coupling shows a logarithmic pole of the form 

e(A) ~ -log(A- A c ). (3.18) 

Hence the mass of the constituent BPS monopoles liberated by dissociating nontrivial holon- 
omy calorons, given as 

47T 

-^monopole ~ ~T , (3.19) 
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Figure 3.1: Evolution of the effective gauge coupling e in the electric phase for SU(2) (grey 
line) and SU(3) (black line). The gauge coupling diverges logarithmically at A c = 11.65 
(SU(2)) and A c = 8.08 (SU(3)). The plateau values are e = 5.1 (SU(2)) and e = 4.2 
(SU(3)). The graph is taken from 0. 

approaches zero with the consequence that magnetic monopoles condense. 

In total, the theory is seen to have three phases: the electric phase at high temperatures, 
the magnetic phase for a small range of temperatures comparable to the scale A, and a 
center phase for low temperatures. The electric phase is deconfining, the magnetic phase is 
preconfining, and the center phases completely confining. Here, we are only addressing the 
physics in the electric phase. 

3.2 The composite adjoint Higgs field <\> 

In this section, the composite adjoint scalar field = a A a (A a denote the Pauli matrices) 
which represents the topologically nontrivial, BPS saturated and trivial-holonomy part of 
the ground state is to be computed. 

The field can be written as a product of modulus and phase. As we are aiming for a 
description of a thermodynamical ground state, and the modulus \(f)\ is a (gauge invariant 
and hence) physical quantity, it has to be homogenous both in space and time. To be 
able to calculate the modulus |0|, the Yang-Mills scale must be known. The phase ^ is a 
dimensionless quantity, so that for its calculation no information about the Yang-Mills scale 
A is necessary. The classical caloron action S = is independent of temperature; this 
excludes explicit /^-dependences of 0's phase. The phase may depend on the temperature 
only via the periodicity of the caloron. Briefly, we can write 



The calorons generating the field are BPS saturated; as an immediate consequence, the 
energy-momentum tensor vanishes identically on a caloron. As already mentioned, the field 
will not include interactions, so the energy-momentum tensor has to vanish on as well. 
This in turn can be described by a BPS equation for 0, 




(3.20) 



«9 T = , 



(3.21) 
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where V^l 2 ^ is the 'square-root' of the potential V = tr^ 1 / 2 ))" 1 " V^ 1 / 2 ). 

In sections 13.2. II and I3.2.3| we will define an adjointly transforming integral over a two- 
point function and demand that 0's phase obeys the same equation of motion. In section 
I3.2.4l the BPS equation for will be determined. In section 13. 2. 51 we will see that introducing 
the Yang-Mills scale A externally allows us to obtain 0's modulus in terms of A and (3. With 
phase and modulus known, the BPS equation is used to determine the potential governing 
the dynamics of 0. The potential is unique, if analytical dependence of the right-hand side 
of the BPS equation on the field is demanded. 

3.2.1 Definition 

We define the phase of the field as proposed in [3] by 

T 







d 6 x / dptiX a F^(r,0) {(r,0),(r,aj)} F^(r,x) {(r, x), (r, 0)} , (3.22) 



where 



tr< 



{(t,x),(t, 0)} = Vexp 



(t,x) 



(r,0) 



(3.23) 



(t,x) 



dZpA^z) 



The Wilson lines in Eq. (13.23)1 are to be calculated along the straight line connecting the 
points (r, 0) and (t,x). V denotes the path-ordering symbol, and tr the SU(2) trace. In 
(I3.22J) . the dependences of the integrand on the caloron scale parameter p and inverse 
temperature (3 suppressed. 

Eq. ()3.22|) is a definition of the phase ^ in the following sense: We demand that ^ 
obeys the same homogenous evolution equation in r as the right hand side of ()3.22|) does, 



V 



101 







(3.24) 



where T> is a differential operator in r. Thus ()3.24j) is an equation of motion for 0's phase. 

The right-hand side of ()3.22|) is to be evaluated both on the caloron-field and the 
anticaloron-field and afterwards the sum is to be taken. We will see in the course of the 
calculation, that the definition ()3.22|) contains quite a number of ambiguities which span 
the solution space of the differential operator T>. 

Under a gauge transformation f2(r, x), the involved objects transform as follows: 



{(r, 0), (r, x)} - fit( r , 0) {(r, 0), (r, x)} fi(r, x) 
{(r, aj), (r, 0)} ^ fit( r , x) {(r, x), (r, 0)} fi(r, 0) 
F Mi/ (r, x) -> fi f (r, x) F^r, x) SI(t, x) . 

Hence the right-hand side of ()3.22|1 indeed transforms like an adjoint scalar, namely 



(3.25) 



h-Xr)^R ab {r)^-{r 



(3.26) 
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where R ab is the SO (3) matrix 



R ab (r)X b = Q(t, 0)A a fi t (r, 0) . (3.27) 

The field <ft only transforms under the time dependent part of the gauge transformation, the 
spatial dependence of the gauge transformation is lost in the macroscopic description. 

The construction in ()3.22|) is invariant under spatial translations and hence the integra- 
tion over spatial translations is trivial. For a gauge variant density as the integrand in ()3.22|) 
is, averaging over global color rotations yields zero and thus is forbidden. The same applies 
to time translations. The only modulus of the caloron solution that is integrated over is the 
scale parameter p with flat measure. 



3.2.2 Are there alternative possibilities for the definition of 0's 
phase? 

The definition ()3.22)1 for 0's phase is not at all arbitrarily chosen. Indeed, trying to generalize 
the right-hand side of ()3.22j) always requires the introduction of either explicit temperature 
dependences (which are not allowed) or additional scales (which on the classical level do not 
exist): 

• Every local definition including the field strength only, such as 

tr A F^ u F UK F Kfl , tr A F^Fj^^F^pFp^ etc. (3.28) 



or 



tr EabcFpvX F UK \ C F KfJ- etc. (3.29) 

yields zero when evaluated on the (anti) caloron field. This is due to the (anti) self- 
duality of the (anti) caloron field Eq. ([2.98)1 . 

One could be tempted to consider higher n-point functions of the type employed in 
Eq. flHI221), such as 

/T 1 / d 3 x j d?y J rfptrA a F^(r,0) {(r, 0), (r, x)} F UK (r,x) 

■{(r,x),(r,y)} F Kfl (r,y) {(r, y), (r, 0)} 

In such an n-point function, one has to introduce an additional factor (3 2 ~ n to get 
a dimensionless object. 0's phase is supposed to depend on temperature only via 
the temperature dependence of the caloron field, but not explicitly. (This is due to 
the temperature- independence of the caloron action.) Therefore the n-point function 
(j3.30J) and its generalization for n > 3 are forbidden. 

One could think of shifting the spatial part of the starting point of the construction 
in ([3.22J1 from to an arbitrary point z. Any given value of \z\ ^ would introduce 
an additional scale; but, on the classical level, such a scale does not exist. The 
same argument applies to replacing the straight Wilson lines with curved arcs. Their 
curvature again corresponds to an additional scale which is physically not present. 
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• The definition (J3.22)) does not include caloron solutions with topological charge \Q\> 
1. This has the following reason: A caloron of topological charge \Q\ > 1 has m > 1 
dimensionful moduli. For example, a caloron with \Q\ = 2 has three dimensionful 
moduli, namely two scale parameters and the distance between its two centers. An 
n-point function of the type displayed in Eq. ()3.30|) contains n field strength tensors 
(mass dimension 2) and n — 1 integrations over 3-space (mass dimension —3). The m 
dimensionful moduli of the caloron have to be integrated as well (mass dimension — 1 
each). These combine to mass dimension 2n — 3(n — 1) — m = 3 — n — m. We are 
looking for a dimensionless object without any explicit dependence on (3. But this is 
not possible with n > 2 and m > 1. Therefore, calorons of higher topological charge 
are excluded. 



3.2.3 Calculation 

In this section, we want to evaluate the right-hand side of 1)3.22)1 . 

J d 3 x j dptr\ a F^(r,0) {(r, 0), (r, *)} F^(r,x) {(r, x), (r, 0)} , (3.31) 
on the single caloron solution (cf. Sec. 2) 

A M (T,aO = -^y^mn(r,r) (3.32) 

with the pre-potential IT(t, r) given as 
vr p 2 sinh 2gE 

n (r, r) = 1 + -^ —c^ " 2nT , (3.33) 



/3r cosh ^ — cos ■ g 



and r = \x\. 



The Wilson lines. The single caloron solution has a hedgehog like behavior in the sense that 
the spatial part of the scalar product x^A^(t, x) has the same orientation both in 3-space 
and color-space, 

Xi A1(t, x) = -x { fjf u d u In U(t, r) 

= -Xi (e aiu - 5 uA 5 ai ) d v In U(t, r) 

= Sain^— 1 d r In Il(r, r) + x { 5 ai 9 4 In II (r, r) 1 " ' ' ' 

r 

= x a (9 4 lnll(r, r) . 



This property allows to discard the path ordering operation in the calculation of the Wilson 
line {(r, 0), (r, x)} defined in Eq. ()3.23j) . We parameterize the path as z^(s) = (r, sx) with 
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< s < 1. Hence we have 

{(r,0),(r, x)}=Vexp 

= V exp 
= "Pexp 



(t,x) 



(r,0) 
1 







exp 



dsxi Ai(r, sx) 



ds ^ Xi <9 4 In II(r, sr) 



ds <9 4 lnn(r, sr) 



(3.35) 



■ x ■ 

cos (g(r, r)) + i A* — sin (g(r, r)) 
r 



where we define 



9( T , r )= I ds-d 4 \nIl(T,sr) . 



(3.36) 



The evaluation of the Wilson line on the single anticaloron solution yields just a change of 
sign in the argument of the exponential function. So, for caloron and anticaloron we have 



{(r,0),(r,x)} 



Xi 



C,A 



cos (g(r, r)) ± i\— sin (g(r, r)) 



(3.37) 



respectively. The following relations hold: 

{(r, 0), (r, *)} c = {(r, x), (r, 0)} A = {(r, 0), (r, -x)}^ = {(r, 0), (r, -*)} 
= {(r, x), (r, 0)} f c = {(r, 0), (r, x)}^ = {(r, -x), (r, 0)} c . 

The integrand of Eq. ()3.36|) for large r behaves like a ^-function in s, 

r 

lim - <9 4 In II(r, sr) = 5{s) ■ f(r) . 

r— >oo 2 

This property has been established numerically. 

Caloron field strength. The field strength F^ u on a caloron can be calculated as 

F; u = d,A a v - d v Al - i [A„ A v ] a (3.40) 

a (d K u)(d K u) a n(d M <9 K n)-2(«9 M n)(«9 K n) _ Q n(9An)-2(p)(a K n) 



(3.38) 



(3.39) 



77 



*7 



+ *7 



The calculation is performed in Appendix IA. 21 To obtain the field strength of an anticaloron, 
77 has to be replaced with t]. In particular, the field strength at x = is 

[5 4 n(r,0)] 2 2<9fn(r,0)' 



^CM) = < 



n 2 (r,0) 3 n(r,o) 



(3.41) 



T/ie integrand of Eq. (|3.31|) . Inserting the result for the Wilson lines Eq. ()3.37|) into the 
expression ()3.3H) and writing the field strength in components, F„ u 



?a 

fit/ 2 ' 



we see that 



trA a F^(r,0) {(r, 0), (r, x)} i?(r, aj) {(r, x), (r, 0)} 



Caloron 



1 

2 tr 



A a A b ( cosg(r, r) + i\ c — sin g(r, r) ) A d ( cos g(r, r) — i\ e — sin (?(r, r 



..r 



,1 



(3.42) 



•F%,{t,0)F*,(t,x 
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is to be computed. Performing the trace and contracting Lorentz and color indices is a 
straightforward but somewhat lengthy calculation; for details see Appendix IA.3I The result 

3 



is' 



trA a F^(r,0) {(r, 0), (r, x)} F^{r, x) {{r, x), (r, 0)} 
.x a /[<9 4 n(r,0)] 2 2 9|n(r,0) 



Caloron 



2i- 



r \ n 2 (r,0) 3 n(r,0) 

2co S (2 9 ( T , r ))f 2 B n < T ' r » n < T -''>l SAn ' T ' r » 



(3.43) 



n 2 (r,r) n(r,r) 

+ sln(2s(T , r)) , 2 lW-f _ M + ■■) r ; 



n 2 (r,r) n 2 (r,r) n(r,r) n(r,r) 



The factor 



[9 4 n(r, 0)] 2 2 «9 4 2 n(r, 0) _ 16vr 4 p 2 *V + P ( 2 + cos 



2ttt 

(9 



n 2 (r,0) 3 II(r,0) 3 /? 2 ( 2 , V + ^ ( x _ cos 



(3.44) 



arises from the field strength at the point (r, 0), cf. Eq. 1)3.41)1 . and contains no dependence 
on r. Note that the expression ()3.43)) is proportional to the (spatial) unit vector 
Integration over position space. The integration over position space demanded in (j3.31j) is 
performed in polar coordinates. On the one hand, the angular integration over the unit 
vector — yields zero. But on the other hand, the radial integral is infinite. This can be seen 
as follows: The pre-potential ()3.33j) for large r behaves as 

II(t, r) - 1 + ^ (r » (3/2tt) , (3.45) 
and its second spatial derivative approaches 

9 2 n(r, r) - dl (l + = (r » /3/2tt) . (3.46) 

Hence the very last term in Eq. (|3.43|) has the asymptotic behavior 

&)^^ (r»m*) (3.47) 
Il(r, r) pr 6 

and gives rise to the logarithmically divergent integral 4 

a 2 nfr,r) .. r°° , . , „2tt 



f°° d 2 Yl{T,r) f°° ! 

/ dr r 2 sin(2g(r, r))— — ^— = finite + / dr sin(2g(r, r))- 
7o n(r,r) y fl /3r 

2vrp 2 / . , ,,\ Z" 00 dr 

= finite H — lim sin(2g(r, r)) / — 

p > Jr r 



(3.48) 



3 Note that d r denotes the derivative with respect to the radial coordinate r = \x\. For more notations 
and conventions, see Appendix lA.il 

4 Note that the integral is convergent at r — 0, since we have 

d r 2 II(T,r) _ 4ttV 2 1 ( 2 + cos (T)) 



r^o n(r, r) 3/3 2 sin 2 27r2p2 cos " 
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where R ^> j3/2n. The function g(r) defined in Eq. (jH.Hfij) has a finite limit for r — > oo, 
which is reached very rapidly. It has been taken out of the integral. 

All the other terms in Eq. ()3.43j) give rise to finite integrals in r, and hence do not 
contribute: For the square of the first spatial derivative, we have 

which is convergent. All the other terms contain at least one time derivative and thus vanish 
exponentially, e. g. 

n , \ 2ixp 2 (2tvt\ ( 2vrr\ „ . N . . 

d T U(r, r) -> —f- sin — -J exp — (r > (3/2n) . (3.50) 



Pr \ P J V P 
Thus the only nonvanishing contribution to the expression (|3.31|) is 



IP(t,0) 3 n(r,0) J Js 2 r J Q v - ' - n (r, r) 

4 7T 2 p 2 +/? 2 (2 + COS^ 1 ) /• o f«)J 

dp^- ) L> rffi / -sin(2,(r,r)). (3.51) 

^ ^ttV + ^^I-cos 2 ^)) ^ r ^i? r 



64i7r 



Regularization. As already pointed out, this expression contains a product zero x infinity 
and thus needs to be given a finite value by prescribing a regularization procedure. The 
radial integration in Eq. (j3.51|) is regularized according to 

dr nP f°° dr , . 

P £ -TT7 (3.52) 



r r J R r x + e 



with e > 0. This integral is 



oo 



The right-hand side is a regular expression for e ^ 0. It can be regarded as the analytical 
continuation of the integral on the left-hand side for |e| <C 1. When smearing the regulator 
e over a small interval [—77, rf\ with < 77 <C 1 as 

1 r . / 1 , {R\ 1 , 2 2 \ 7T7 /iA 2 



2, * (7±Io - los UJ + ? los U ) + ) = T ^ " '° g U J + a " ' (3 ' 54) 

an ambiguity appears: There are two possibilities how to circumvent the pole. 

The angular integration is regularized via the introduction of a defect (or surplus) angle 
in the azimuthal integration, 

ddsmO / dif^ d6 sin 6 dtp (3.55) 

JO JO Ja c ±v' 
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with < rf <C 1 and an arbitrary angle < ac < 2ir. The regulated radial integral yields 
I d8 sin 6 dip — = ^fn smr]'(5 al cosa c + 5 a2 sma c ) , pr^?\ 

JO Ja c ±v' T (3.56) 

= =R"/(£»i cos «c + ^2 sin a c ) + 0r]' 2 ■ 

The regularization prescription Eq. (J3.55)) clearly singles out the axis with unit vector 
(cos etc, sin etc, 0). Later this will be seen to relate to a global gauge choice. 

At the end of the calculation, both regularization parameters r\ and rf have to be sent 
to zero. Their ratio in this limit will be denoted by Sc, 

lim - = E C . (3.57) 

This is a positive but otherwise unknown number. Inserting the regularized expressions 
Eqs. fling) and ()3~5fij) and into ((HEED , we get 



J d 3 x J dpti\ a F^(T,0) {(r,0),(r,aj)} F^{r,x) {{r, x), (r, 0)} 
= ±H c (5ai cos a c + S a2 sin ac)^4 



Caloron 



2vrr 



(3.58) 



/3 

where the dimensionless function A is defined as 

A f2nr\ 327T 7 /• / , A 4 vr V + /3 2 (2 + cos ^ , 
.A — - = — — / dp lim sin2^(r,r) p 4 - ^ f-^ . (3.59) 



(3 J 3 /? 3 J r V™ 



(2vr 2 p 2 + (3 2 (l - cos 



2ttt 

& 



Integration over scale parameter. Up to now, we did not specify the range for the integration 
over the scale parameter p. As the integrand in Eq. ()3.59|) asymptotically behaves like ~ p 2 , 
integrating p from zero to infinity yields an infinite expression. To see what is going on, a 
cutoff is introduced in units of inverse temperature /3, 

rCfl 

dp^ dp (C > 0) . (3.60) 



This generates an additional dependence of A on (. For £ ^> 1, we will have A oc ( 3 . The 
integral in Eq. ()3.59|) has to be evaluated numerically. Due to the property Eq. ()3.39|) . the 
limit lim r _ i . 0O sin 2g(r, r) is reached very fast; for our purposes, putting r = 10 in numerical 
calculations is fully sufficient. Fig. 13.21 shows the r-dependence of A for various values of £. 
With growing (, the function A rapidly approaches a sine curve, 

a(^x) ^ 272C 3 sinf^) (C^oo). (3.61) 



P 'V * \ (3 

Already for ( = 10, the difference between the two curves could not be resolved any more 
in the figure. The prefactor 272 has been fitted numerically, see Table 13.11 Therefore, the 
evaluation of the expression 1)3.31)1 on the single caloron yields 



d 6 x J dp trA a F^(r,0) {(r, 0), (r, x)} F^(r,x) {(r, x), (r, 0)} 
±272C 3 S c (5 al cos a c + S a2 sin a c ) sin {—- 



Caloron 



(3.62) 
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Figure 3.2: The function A(-j L ) plotted over two periods with different values of £. For 



comparison the function 272( • j sin(^) is plotted dashed line. Already for £ 
difference cannot be resolved any more. 
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Table 3.1: Val 

dependence ( 3 has been divided out. 



c 




C 


A(v/e 


c 


A(v/e 


1 


301.295 


10 


272.776 


100 


272.026 


2 


285.012 


20 


272.216 


200 


272.020 


3 


278.828 


30 


272.107 


300 


272.018 


4 


276.161 


40 


272.068 


400 


272.018 


5 


274.794 


50 


272.050 


500 


272.018 










1000 


272.018 


of the function A at 


| for several values 


of the cutoff 



3.2.4 Ambiguities and BPS saturation 

The expression f!3.31|) has been evaluated on the caloron in the previous section. Evaluation 
on the anticaloron yields the same result, if we agree upon circumventing the pole appearing 
in Eq. ()3.54j) in the opposite way as compared to the caloron. According to the definition, 
both contributions are to be added. In this process some ambiguities occur: 

1. The undefined number lim^ ir) '-»o — needs not necessarily be the same in the caloron 
and anticaloron case, we chose to call it for the anticaloron (instead of for the 
caloron) . 

2. The same applies to the arbitrarily chosen axis singled out by angular regularization. 
It is no restriction of generality to assume that both axes lie in the xia^-plane, but 
with different azimuthal angles «c and «a- 

3. As ()3.22|) is supposed to define only the equation of motion for 0, we may as well 
introduce a shift in time; this does not change the operator T>. Again, this shift 
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may be different for the two contributions, namely r — > r + Tc for the caloron and 
t — > r + ta for the anticaloron. 

The requested sum of caloron and anticaloron contribution (including all the above ambi- 
guities) is then 



— ~ 272<" 3 <( ±S c (S a icosa c + S a2 sma c )sm 
± E A (6 a i cos a a + S a2 sin a a) sin 



2tt(t + T C ) 

(3.63) 

2n{T + T A y 

(3 



The caloron (or anticaloron) contribution in Eq. ()3.63|) alone represents a linearly polarized 
harmonic oscillation in adjoint color space. The three ambiguities given above can be viewed 
as the free parameters of such an oscillation, namely modulus, phase-shift and polarization 
axis. The sum of caloron and anticaloron contribution Eq. (J3.63)) taking into account all the 
above ambiguities is an elliptically polarized oscillation in adjoint color space. Here, the po- 
larization plane is the XiX2-plane; this is only due to our choice of the angular regularization 
and no physical property. 

The right-hand side of Eq. ()3.63|) obviously is annihilated by the second order differential 
operator 

V = dl+( 2 4\ ■ ( 3 - 64 ) 



P 

Note that the ambiguities inherent to the definition ()3.22j) span the solution space of T>. 
Imposing BPS saturation. The field (f) is composed of noninteracting trivial-holonomy 
calorons; the caloron itself is an energy- and pressure-free, BPS saturated configuration 
and no interaction whatsoever has been included. Hence the composite object must again 
be energy- and pressure-free and BPS saturated; that means it does not only obey a second 
order differential equation (the equation of motion), but also a first order differential equa- 
tion (the BPS equation). Thus we need to find first-order equations whose solutions solve 
the second order equation 

afr+(j)% = (3.65) 

as well. There are two such equations 5 , namely 

d T = ±^A 3 0. (3.66) 

Choosing any (normalized) linear combination of Pauli matrices instead of A3 would be 
possible equally well: Eq. ()3.66j) is subject to a global gauge ambiguity. The solutions to 
()3.66|) are given as 

= CAi exp (±^A 3 (r - 7-0)) , (3.67) 



5 The solutions to the equations d T (f> = ±^(j) as well solve Eq. H3.65[) : in spite of this they are not allowed 
here because they are not adjoint fields. 
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where C and To are real integration constants. This solution is a circularly polarized oscil- 
lation. It winds along an S\ in the group manifold S3 of SU(2). 

Thus the demand for BPS saturation forces an elliptical polarization in Eq. ()3.63j) into 
a circular polarization. The undetermined and formerly independent quantities are now 
subject to the relations 

71 7T 

■=c = -a , ro - Ta = ±- , ol c - ol a = ±- . (3.68) 

The modulus of the oscillation and its phase-shift are still undetermined constants of inte- 
gration. The former will be considered in the following section, the latter is of no physical 
significance. In addition, a global gauge ambiguity, i. e. the plane in which the oscillation 
takes place, is still present. 



3.2.5 Obtaining 0's modulus 

Let us now assume the existence of an externally given scale A which determines 0's modulus. 
We allow for explicit dependence of on the scale A, the temperature j3 and on Euclidean 
time through 5, 

= (p, A, lj . (3.69) 

As the phase found in Eq. ()3.67|) shall be preserved even in case of the presence of a scale 
A, the right-hand side of the BPS equation 

d r( f) = V {1/2) (3.70) 

may only depend linearly on 0. Besides that, we demand an analytical dependence of V on 
0. The potential V (and its 'square- root' V {l/2) ) may depend on the temperature through 
the periodicity of 0. An explicit dependence on (3 is not possible because no explicit (3- 
dependence occurs in the average over the caloron moduli space. 
These conditions leave only the two possibilities 

<9 T = ±«AA 3 (3.71) 

and 6 

«9 T = ±zA 3 A 3 A. (3.72) 
Using Eqs. ()3.2()j) and ()3.67|) . we write 

0= |0(/3,A)|A ieX p (±^pv) • (3.73) 

The modulus of is gauge invariant and hence a physical quantity describing a thermo- 
dynamical ground state, so it has to be homogenous in space and time. Inserting the 
decomposition ()3.73|) into the first BPS equation (|3.7ip . we get 

A = j ■ (3-74) 



Note that j^p = (j> 1 = O 1 S^Lo(~ l) n 0o "(^ — ^o) n is indeed an analytical function of 1 
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This obviously can not be satisfied, since the scale A is a constant and (3 is the inverse 
temperature. From the second possibility Eq. ()3.72j) . we get 



which is no contradiction. So, Eq. (|3.72j) is the only acceptable BPS equation, and has 
the form 



/3A 3 [ 2m 

Axexp ( ±— — A 3 r . (3.76) 



2?r r V P 

Eq. ()3.75|) shows that the field is power suppressed in T, and hence all topologically 
nontrivial effects die off at high temperature. The right-hand side of the BPS equation 
defines the 'square-root' of 0's potential, 

V{$) = tl (V^yV^ = A 6 tr0~ 2 . (3.77) 

Under the above assumptions, the potential is unique. The Lagrangian for the field is 

C = tr(<9 T 0) 2 + 1/(0) . (3.78) 

The consequences of minimally coupling to the topologically trivial sector have been 
investigated in [H]; some of the results are briefly reviewed in Sec. 3.1. 
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Chapter 4 



Thermodynamical Pressure in the 
Electric Phase 



In this chapter the two-loop corrections to the thermodynamical pressure in the electric 
phase of SU(2) Yang-Mills theory are computed. In view of the evolution of the effective 
gauge coupling in the electric phase (cf. Fig. l3.1|) . one can constrain oneself to the case e > ~. 
This will simplify the calculation. For a numerical evaluation, the plateau value e = 5.1 
is used. We will set up the prerequisites for the calculation in Sec. 14.11 The calculation is 
performed in Sec. 14.21 In Sec. 14 .3\ the two-loop contributions are compared to the one-loop 
result for the pressure. 



4.1 Feynman rules and other prerequisites 
4.1.1 One- loop pressure 

The thermodynamical pressure P is defined as the derivative of the partition function Z 
with respect to the volume of the system, 

In a field theory, P can be calculated order by order in a loop expansion with the help of 
diagrammatic techniques as presented in [2*B1 12bj . 

The pressure has been calculated on the one-loop level for SU(2) and SU(3) Yang-Mills 
theory in [Sj; we only consider the SU(2) case here. The pressure contains a temperature 
dependent ground-state contribution arising from caloron 'condensation', 

P g . s .(A) = -2AA 4 , (4.2) 

and a contribution associated with the one-loop diagrams shown in Fig. 14.11 

Q\4 / 27T 4 - \ 

Pone-loop(A) = -A 4 ■ j^y 6 + 6P ^A^J . (4.3) 

The fixed mass scale A is connected to the modulus of the composite Higgs field \<f>\ by 

W 2 = ^, (4.4) 
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Figure 4.1: One- loop contributions to the pressure. Double lines correspond to tree- level 
heavy (TLH) modes, and single lines correspond to the tree-level massless (TLM) mode. 



the dimensionless temperature A is defined as 
2vrT 



A 



A ' 

and the function P is given as 



P{a) 



dx x 2 log 







1 — exp ( — \Jx 2 + a 2 



(4.5) 



(4.6) 



There is also a 'nonthermal' contribution —AV. It is estimated as 
AV < 



16tt 2 



Because 

AV 



V 



< 



32vr s 



< 2 • 10" 



(4.7) 



it can be neglected in the electric phase. (Recall that A > 11.65 in the electric phase, cf. 
Sec. 3.1.) 




4.1.2 Feynman rules 

For the calculation of the two- loop correction AP to the pressure of SU(2) being in its 
electric phase, we have the equation (see jH]) 



(4.9) 



where double lines represent TLH modes and single lines stand for TLM modes. TLM 
modes will carry a color index 3, while the color indices 1 and 2 correspond to TLH modes. 
Eq. (j4.9j) is valid for SU(N). In the case N = 2 considered here, the second and third diagram 
do not occur. 

The calculation of the thermodynamical pressure is performed in unitary gauge, where 
(ft is diagonal and the background is a b ® = 0. The remaining gauge freedom is used to gauge 
the TLM mode to transversality, diSaf LM = (Coulomb gauge). Unitary- Coulomb gauge 
is a completely fixed gauge, thus no Faddeev-Popov determinant needs to be considered and 
no ghost fields need to be introduced. 

The computation is performed in the real-time formulation of thermal field theory. 
Spacetime is Minkowskian with signature (+,—,—,—). Matsubara sums, which origi- 
nate when working in compactified imaginary time, are replaced with integrals over real 
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Minkowskian momenta by means of contour integrals and analytic continuation. For our 
purposes, the real-time formalism is preferable because the implementation of constraints 
for the momenta of the participating fluctuations (Eqs. (|4.15jl . (|4.1fi|) ) is rather inconve- 
nient in the imaginary-time formalism. Moreover, the quantum and thermal part of a given 
fluctuation can be clearly discerned in the real-time formalism. 

In the real-time formulation and in unitary- Coulomb gauge, the Feynman rules employed 
in calculating the diagrams in Eq. ()4.9j) are as follows: The propagator for a free TLM-mode 

0E51EE] 

D T ^(k,T) = -5 ab P^(k) (^-— + 2n5{k 2 )n B {\h\/T)^ , (4.10) 
where 1 

pJW = f£M = Pat*) = o , 

pm=s»-f. (4 ' u) 

and n B denotes the Bose-Einstein-distribution, n B (x) = -^zj- ^he P ro P a g a tor for TLH- 
modes reads 

DlZ^T) = -5 ab [g, v - ^ ( k 2_^ + t£ + 2 ^ - - 2 ) n B (\ko\/T^ . (4.12) 

The vertices are the usual ones (see j2HJElj): The four-boson- vertex is 
fi, a v, b 

XX = -*e 2 (2vr) 4 5\p + q + r + s) \e jab e Scd ^g v ° - g^g^) 

S /\ +e fac e fdb (g^g^ - g» v g»°) + e fad e fbc (g^g<"> - g^g™)] , (4.13) 

and the three-boson-vertex is 
v,b 

Q =e(2ir) 4 6\ P + q + k)e abc {g^(q- P y + g^(k-qr + g^(p-ky]. (4.14) 




According to Landsman and van Weert |2Sj , one has to divide every diagram by % and by the 
number of its vertices. Besides that, we demand that the momenta of quantum fluctuations 
are to be cut off at 

\p 2 ~m 2 \ < |0| 2 (4.15a) 
in Minkowskian or 

p% + m 2 < |0| 2 (4.15b) 



1 Static electric fields of long wavelength are completely screened due to the existence of an infinite real 
part in the Debye screening mass mo = [IIoo(&o = 0, k — > 0)] 1 / 2 ; cf. App. B 
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in Euclidean signature because of the existence of the compositeness scale \<j>\, cf. Sec. 3.1. 
The compositeness scale imposes a similar constraint on the center-of-mass energy flowing 
into the four- vertex, 

|(pi+p 2 ) 2 |<H 2 , (4.16) 
where pi are the ingoing momenta. 



4.1.3 Contributing diagrams 

The various contributions to the pressure on the two-loop level will be denoted as follows: 
The first (last) diagram in Eq. ()4.9|) shall be referred to as AP MHH (AP HH ), since there 
are two TLH and one TLM particle (two TLH particles) involved. If we want to specify, 
for example, that in AP MH the TLH fluctuation is a vacuum fluctuation and the TLM 
fluctuation a thermal fluctuation, we write AP^y 1 etc. The statistical factors of Eq. (|4.9|) 
are included in these definitions. 

Which diagrams in Eq. (|4.9jl do actually contribute for SU(2)? Because of the structure 
of the vertex f)4.14j) , the second and the third diagram in Eq. (0~HJ vanish. Moreover, AP^™ 
vanishes (for any N and e), since the on-shell conditions p 2 = k 2 = m 2 , q 2 = and energy- 
momentum-conservation p + k + q = cannot be satisfied simultaneously. 

As the mass of TLH fluctuations is connected to the gauge coupling via m 2 = 4e 2 \4>\ 2 , 
the constraint (|4.15jl reads 

|p 2 -4e 2 |0| 2 | < |0| 2 or p\< (l-4e 2 )|0| 2 . (4.17) 

This condition will eliminate all contributions containing quantum fluctuations of TLH 
modes, if the gauge coupling e is larger than 1/2. In this case, we only have the following 
nonvanishing contributions: 

AP HH = AP™ 

AP MH = AP^ T H + APyJjP (4.18) 
AP MHH = AP™. 

So, using the above notation, Eq. ()4.9j) for e > 1/2 simplifies to 

AP = AP™ + AP$? + AP^ T H + AP™ . (4.19) 



4.2 Calculation 

4.2.1 Calculation of AP™ 

We first write down the expression for AP™ without imposing the kinematical constraints 
Eqs. (J4.15|) . ()4.16|) and perform the constrained integrations afterwards still using the nota- 
tion AP™. According to the Feynman-rules given in section l4.1.2[ the diagram shown in 
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H, a p, c 



ap hh = _.p 

8 




v . b a, d 



Figure 4.2: A local bubble diagram containing two TLH-modes 



Fig. 14.21 with both fluctuations being thermal reads as: 



+e abf e fdc (g^g»» - g^gT) + e acf e fdb (g^g^ - g^g'")] 

■ (Sab) (g»» - 2n5(k 2 -m 2 )n B (\k \/T) 

■ (-**) (g P „ - P -if) 2vr 5{p 2 - m 2 ) n B (\p \/T) 

p 2 k 2 p 2 k 2 (pkf 



d 4 k f d 4 p 



12-3^-3 — + 



m 2 m 2 m 4 m 4 



2 J (2vr) 4 J (2tt) 4 
2tt 5(k 2 - m 2 ) n B (\k \/T) ■ 2ir5(p 2 - m 2 ) n B (\p \/T) . 



(4.20) 



(The contraction of Lorentz and color indices is deferred to Appendix IB. 11 ) In evaluating 
this integral, we will first perform the integrations over the zero-components of the momenta 
and thus eliminate the 5-functions. Therefore we rewrite the 5-functions as follows: 



AR 



HH 
TT 



e 2 f d 3 k f d 3 p 



2 J (2tt) 3 J (2tt) 3 
1 



dk / dp [7 



p\k 2 - 2p k pk + (pky 



2Vk 2 + m 2 



5 (k - Vk 2 + m 2 J +8 (k + V k 2 + m 2 
S (po - a/p 2 + m 2 ) + 5 (po + \/p 2 + rr 



■n B {\k \/T) (4.21) 
• n B (\p \/T) . 



Now we notice that the two terms containing 



5 (k — V k 2 + m 2 ^ 5 (p — \/ p 2 + m? 



(4.22) 



or 



5 (k + Vk 2 + m 2 ) 5 (p + vV + m 2 ) 



(4.23) 



respectively, yield the same result since the rest of the integrand is invariant under simul- 
taneous reflections po — ► —po, k — > — k . The same holds true for the remaining two 
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summands. So the expression is split into two contributions as 

ap hh = -—^r[ d3k [ d3p 1 

p 2 + m 2 ){k 2 + m 2 ) =F 2a/p 2 + m 2 \/k 2 + m 2 pfc + (pfc)' 



7- 



m 4 



n B (Vp 2 + m 2 /r) • n B (V& 2 + m 2 /T) 
2 r d 3 k f d 3 p 1 



e 



(4.24) 



4 (2^)3 7 ( 27r )3 v / fc 2 + m 2 v / p2+m2 

p 2 fc 2 p 2 fc 2 + (pfc) 2 a/p 2 + m 2 a/A; 2 + m 2 
6 ^ ^ 4 ± 2pfc r 

m z m z m 4 m 4 

• n B (^p 2 + m 2 /T^ ■ n B (^k 2 + m 2 /T^j . 

The tree-level mass of TLH-modes in SU(2) is given as 

m = 2e\<j>\, (4.25) 
so the introduction of dimensionless momenta y and x according to 

k = y ■ \4>\ and p = x ■ \<f>\ (4.26) 

is suggested. Moreover, temperature T and compositeness scale |0| can be expressed in 
terms of the dimensionless temperature A and the mass scale A (see Eqs. (j4.4j) and (|4.5jl ) as 

|0| = AA~ 1/2 , 

(4 27) 

|0|/T = 2ttA- 3/2 . 

Finally, we introduce three-dimensional polar coordinates for the scaled momenta. The 
angle between x and y can be chosen to be the polar angle 6. AP™ is then expressed as 

e 2 A 4 A- 2 ^ r , r , f 1 , „ x 2 y 2 



e A A /"°° Z" 00 /" 



± 

2 



\Ac 2 + 4e 2 yV + 4e 2 



I/ 2 h . 2m , o a Vx 2 + 4e 2 y/y 2 + 4 2 \ (4.28) 

— 7 1 + cos ± 2xy cos 6 1 -. v ' 

4e 2 4e 2 16e 4 v 1 y 16e 4 J 

■ n B (2tc\- 3 / 2 Vx 2 + 4e 2 ) • n B (27rA" 3/2 vV + 4e 2 ) . 

Now we have to implement the constraint for the center-of-mass energy in the vertex, 

\(p + k) 2 \ < \(j)\ 2 . (4.29) 



Of course, this condition has to undergo the same manipulations as the integrand itself: p 



2 



and k 2 have to be replaced with m 2 , pok with ±a/p 2 + m 2 Vk 2 + 



m 2 



2m 2 ± 2^p 2 + m 2 Vk 2 + m 2 - 2pk < |0| 2 , (4.30) 
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the momenta have to be rescaled as Eq. (|4.2fi|) . the expression for the mass from Eq. (|4.25j) 
has to be inserted, and a change to polar coordinates has to be done. We then have 

4e 2 ± Vx 2 + 4e 2 vV + 4e 2 - xy cos 9 <-. (4.31) 

For e larger than condition (|4.31|) with the positive sign cannot be fulfilled, and the 
corresponding integral in Eq. (|4.28J) is zero. For the negative sign, condition ()4.31|) can be 
fulfilled, and we will implement it as an additional constraint on the region of integration 
for cos#. The condition 

4e 2 - Vx 2 + 4e 2 vV + 4e 2 - xycos9 = -4e 2 + Vx 2 + 4e 2 vV + 4e 2 +x?/ cos6» < - (4.32) 

is equivalent to 

cos#< — ( - + 4e 2 - Vx 2 + 4e 2 vA/ 2 + 4e 2 ] = g(x,y) , (4.33) 
xy \2 J 

and thus cos6> is to be integrated in the range [—1,1] H (—00,(7]. This can be done by 
replacing the upper limit of integration for cos 9 with 

max {— 1, min {1, g}} . (4.34) 

So, we have in total 

max{ — l,min{l,g}} 

e 2 A 4 A~ 2 r°° , r°° J r , „ x 2 y 2 

AP$? = — / dx dy / dcos9 



2(2vr) 4 Jo Jo J Vx 2 + 4e 2 ^y 2 + 4e 2 

x 2 y 2 xV v^^VF+^V 4 - 3 ^ 
^-^-^- T ^( 1 + cos ^)-2xycos^ — 

■ n B (27rA- 3 / 2 v / a; 2 + 4e 2 ) • n B ^ttA" 3 / 2 vV + 4e 2 

In this form the integrals can be performed numerically. Instead of doing this, we will have 
a closer look at the function g in order to be able to give the boundaries of integration 
explicitly. For the following compare with Fig. 14. 31 

If g is greater than 1, Eq. ()4.33|) is fulfilled automatically and cos# has to be integrated 
from -1 to 1; this is the case in the white area in Fig. 14. 31 bounded by the curve 



-x - 8e 2 x + y/l + 16eV^ 2 + 4e 2 
This curve intersects with the at 



b(x) = — . (4.36) 



1 + < 4 ' 37 > 

If — 1 < g < 1, the upper limit for the integration of cos# is g. The region where this is true 
is shaded in grey in Fig. l4.3j it is enclosed by the lines 



, N x + 8e 2 x -y/l + 16eVa: 2 + 4e 2 

, (4.38) 

„ . x + 8e 2 x + y/l + 16e 2 Vx 2 + 4e 2 

d{x) = 8^ ' 
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y 




Figure 4.3: For determination of integration limits for cos# depending on x and y. Left- 
hand side: In the white region, the upper limit of integration is 1; in the grey region, it is 
g{x,y) as defined in Eq. Ij4.33j) . From the black region no contribution arises. Right-hand 
side: The integration is split into two bounded parts (I and II) and one unbounded part 
(III) as indicated. 

and b(x) as given above. For all other values of x and y (shaded in black), we have g < — 1 
and thus no contribution. This can be summarized as follows: In Eq. (14.35(1 the integral 
f dx J dy f dcos9 is decomposed as 

ra r-b /*1 pa pd pg poo pd pg 

I dx j dy I dcos9+ / dx dy dcos9+ / dx dy dcos8 (4.39) 

</0 JO J-l JO Jb J-l Ja Jc J-l 

where integrand and prefactor are as in Eq. (|4.35jl . The integrations over cos 9 can be 
performed analytically; further evaluation needs to be done numerically. For the result see 
Fig.Pl 

Fig. 14.31 illustrates that the integration over loop momenta is strongly restricted to a 
narrow band around \p\ = \k\ due to the compositeness constraint Eq. (|4.16(1 . 



4.2.2 Calculation of AP^ H and AR 



/i, a p, c 



^pMH 




v . h a, d 



Figure 4.4: A local bubble diagram containing one TLH-mode and one TLM-mode 
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For the diagram depicted in Fig. l4.4[ the Feynman rules yield 



AP 



MH 



d A k 



d 4 p 



-e 2 ) [e ad fe fbc (g^g^-g^g L 



(2tt) 4 J (2tt) 4 
e abf e fdc (g^g up - g w g va ) + e acf e fdb (g^g^ - g^g")] 



- , ; PpPu 

-Ocd) [9 P a - — y 



p2 _ m 2 _|_ l e 



+ 27r6(p 2 -m 2 )n B (\ Po \/T) 



(4.40) 



d 4 k 



k 2 + is 
d A p 



+ 2n5(k 2 )n B (\k \/T) 



(2vr) 4 J (2tt) 



-0 -4 • 4 



i p 2 (pfe) 5 



m 2 m 2 m 2 fc 2 



p 2 — m 2 + ie 



+ 27r5(p 2 -m 2 )n B (|p |/T) 



/c 2 + 



+ 2n5(k 2 )n B (\k \/T) 



Calculation of AP^ T H 

For both fluctuations being thermal, we have 

- 2 <• d 4 k r d 4 P 



APr 



MH 
TT 



r , P P (P*0 

-b + 2 — - H - 



m 2 ' m 2 m 2 k 2 



2 J (2tt) 4 J (2tt) 4 
27t5(A; 2 ) 7i B (M/T) ■ 2n5(p 2 - m 2 ) n B (\p \/T) . 



(4.41) 



The evaluation of this diagram is similar to the calculation in Sec. 14.2.11 Integrating over 
the zero components of p and k, we have 



d 3 k 



2 J (2tt) 3 J (2tt 



d3p f dk J dpo (-4 + 2^ 2 -^ 2 



2 k 2 



2 v P 2 + wz 2 
1 



5 (po - VP 2 + ™ 2 ) + 5 [Po + VP 2 + m 2 ) -n B (\p \/T) 



2Vk 2 

,2 



6 [kn — V k 



6 [k 



n B (\k \/T) 



(4.42) 



e 
~2 



(2vr) 3 J (2vr) 



-4 + 2 



p 2 (pky 



2 fc 2 



2^p 2 + m 2 v / fc 2 



Inserting the mass given in Eq. ()4.25|) . scaling momenta as in Eq. (|4.26|) . introducing dimen- 
sionless variables as in Eq. (|4.27jl and polar coordinates, we obtain 



2 A 4 A- 



2(2tt)- — jo 



/ cfx / cfy / dcos9 



-4 + 



2 a; 2 cos 2 6 



2 



(4.43) 



4e 2 4e 2 J ^ x 2 + 4e 2 
■ n B ^TrA-^Va; 2 + 4e 2 ) • n B (2n\- 3 / 2 y) . 
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At this stage the ± contributions in this equation are the same, but after imposing the 
condition 

\(p + k) 2 \ < \cf>\ 2 , (4.44) 
this is no longer the case. By the above manipulations, Eq. ()4.44|) is transformed into 

±2vV + m 2 Vtf -2pk < |0| 2 (4.45) 
and finally into 

4e 2 ± 2yVx 2 + 4e 2 - 2xy cos 9 <1. (4.46) 

Again, condition 1)4.46)1 . which limits the range of integration for cos 9, can only be fulfilled 
for the negative sign. The condition 



m 2 



4e 2 — 2yVx 2 + 4e 2 — 2xy cos 9 
is satisfied if and only if 



< 1 (4.47) 



-1 + 4e 2 - 2yVx 2 + 4e 2 „ +1 + 4e 2 - 2y^x 2 + 4e 2 , Aln . 

= r 2 < cos 6* < Ti = , (4.48) 

2xy 2xy 

and thus cos 9 is to be integrated over the interval [—1, 1] H [r 2 , Ti). So the upper and lower 
limits of integration in Eq. (J4.4H)) are 

max { — 1, min {1, Tx}} and min {1, max { — 1, T2}} (4.49) 

respectively. Hence the result is 

max{— l,min{l,Ti}} 

min{l,max{~l,T 2 }} (4.0UJ 



X 2 X 2 COS 2 9 

-4 H 

4 e 2 4 e 2 



^ • n B (2n\- 3/2 Vx 2 + 4e 2 ) ■ n B {2^\- z/2 y) 



Again, we want to give the limits of integration explicitly. To do this, the points with 
T1/2 = ±1 need to be determined. They are 

4e 2 + 1 

n = ±1 y = yi± 



r 2 = ±1 <^> y = y 2 ± 



2 (±x + Vx 2 + 4e 2 ) 
4e 2 - 1 



2 , ( 4 - 51 ) 



2 (±x + Vx 2 + 4e 2 ) ' 
The following inequalities are satisfied: 

2e 



V» Vi- « ^" 5 vSF 7T (4 . 52) 

2/2+ < 2/i+ < 2/2- < 2/1- ^ x > r/ . 
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Taking into account all possible combinations of x and y and the corresponding limits of 
integration, the integral in Eq. (|4.5()jl splits into six parts, namely 

v rV2- rl rv rvi+ /•! rv rvi- r T i 

dx dy dt + dx dy dt + dx dy dt 

J y 2 + Jt2 jo Jy2- J—i Jo Jyi+ J—i i a t-q\ 

00 ryi+ r°° fV2- rn roc pyi- pn ^ ' ' 

dx dy dt + dx dy dt + I dx dy dt. 



J V2+ Jt 2 Jr) Jyi+ Jt 2 Jij Jy2- J—l 

The result of the numerical evaluation is shown in Fig. 14.71 



Calculation of APlf H 



VT 



For this case, we have from Eq. ()4.40|) the expression 



A pMH 
^ A - r VT 



e 

T 



d A k 



d A p 



(2vr) 4 J (2tt) 



2n6(p 2 -m 2 ) n B (\p \/T) 



-6 + 2^ + 



p 2 (pky 



m 2 ' m 2 m 2 k 2 J k 2 + is 



(4.54) 



As there is one thermal fluctuation (p) and one vacuum fluctuation (k) involved, only po 
can be integrated by eliminating the 5-function, 



a pMH 
ar VT 



e 

T 



d 4 k f d 3 p 



dpon B (\/p 2 + m 2 /T 



k 2 + is 



p 2 {pk) 2 \ S (Po ~ vV + ™ 2 ) +$(po + \/p 2 + m< 



2 ^ 

± 

-4 + 



m 2 m 2 k 2 
d A k 



2 v P 2 + m 2 
d 3 p ( /— ^ ,„\ % 



(4.55) 



(2tt) 4 J (2vr) 3 
p 2 (pk) 2 \ I 



n B (a/p 2 + m2 /T 



k 2 + is 



i 2 k 2 J 2a/p 2 + m 2 
The implementation of the compositeness constraints 
|(p + A;) 2 |<|0| 2 



and 



k 2 < 



(4.56) 



(4.57) 



is more difficult as in the previous calculations. We will therefore ignore Eq. (|4.5fij) and give 
an estimate for APyJp by only taking into account (|4.57jl . The two contributions ± are 
equal, k is analytically continued to Euclidean momenta. This yields the upper bound 



MHI 
VT 



< 



d 4 k E f d 3 P n B (\/V + m V T ) 1 



(2tt) 4 7 (2tt)3 ^ p 2 + m 2 k 

45 



-4 + 



p 2 (pky 



2 k 2 



(4.58) 



AP 



MHH 




v, b 



«, 9 



Figure 4.5: A nonlocal diagram containing two TLH and one TLM mode 



As usual (inserting the expression for the mass, scaling the momenta and polar coordinates 
for the 3- vector p and the 4- vector fog), we obtain 



2a4\-2 poo 



e 2 A 4 A 



(l.r I dcosti 



X 



K2vr) 4 J J_ x Vx 2 + 4e 2 



4+ 4^ (1 ~ COs2 ' 



• n B (27rA- 3/2 Vx 2 + 4e 2 ) • 
The numerical evaluation of this upper bound is shown in Fig. 



(4.59) 



4.2.3 Calculation of AP^™ 

The contribution AP^jf 1 shown in Fig. 14. 51 corresponds to the expression 
1 e 2 f d 4 k f d 4 p f d 4 q 



AP 



MHH 
VTT 



4 2 % J (2k) 4 J (2k) 4 J (2k)* ^ *<P + « + k ) 
e acf [g^ip - k) x + g<> x {k - qY + g Xfi (q - p) p ] 
e bdg [g™(k - pY + sT{q - k)» + g™(p - qY] 

■ (-cU) (g,u - Vj ^f) 2n5(p 2 -m 2 )n B (\p \/T) 
vpk a 



-Sad \ g P e 

2 r d 4 p 



rn 

knkn\_ „, 2 2\ t\U\lrn\ I X \ tdT 1 



16 I m 



m 

4 - ' d 4 k 



2 2n5(k 2 - m 2 ) n B (\k \/T) ■ (-5 fg ) P^(q) — . 



(4.60) 



(2tt) 4 J (2tt) 4 (k + pY + ie 
j (kp) 2 \ k 2 p 2 -(kp) 2 



77T 



+ 4 



(fcp) s 



(p + kY 

■ 2n 5 (p 2 - m 2 ) n B (\po\/T) ■ 2k 5 (k 2 - m 2 ) n B (\k Q \ /T) . 

The vacuum propagator [(k + pY + is} -1 has a singularity at (p + A;) 2 = 2(m 2 + kp) = 0, 
that is at kp = —m 2 . In spite of this, the term containing 



16 (m 2 



(kpf 



is regular everywhere even for e = because 

16 m 4 — (A;p) 2 16 m 4 — {kp) 2 8 / ,, 



(k + pf 



m 2 2(m 2 + kp) m 2 
46 



— -(m — kp) 



1 _ tE. 



(4.61) 



(4.62) 



So is can be dropped for this term. This is, however, not the case for the second part which 
is proportional to 



fcV - (kpf ( + 1 (fcp) 2 



(4.63) 



In this diagram we again have to take care of constraining the momenta of vacuum fluctu- 
ations, 



\q 2 \ = \(p + kf\ = 2 \m 2 + pk\ < 



(4.64) 



This is exactly the same condition as has been found in the calculation of AP™ (compare 
Eq. ()4.29j) and the following), so we can take the integration limits from Sec. 14.2.11 

We once again reformulate the 5-f unctions, 



a pMHH 
AA - r VTT 



d 3 p f d 3 k 



4 J (2tt) 3 J (2?r) 3 
1 



dp / dk 



2(m 2 + kp) + is 



2 a/p 2 + m 2 



8 [Po- \ / p 2 + m 2 ) +5[po + y / p 2 + m 2 ) n B (\p \/T) 



2Vfc 2 + m 2 
16 ( m 2 



5 (k - Vk 2 + m 2 ) +5 (k + V k 2 +m 2 



n B (\ko\ /T) 



(4.65) 



(kpf\ k 2 p 2 - (kp) 



{p + ky 



+ 4 



(kpf 



and perform the integrations over p and k Q , 



A pMHH 
VTT 



d 3 p f d 3 k 



16 J (2tt) 3 J (2tt) 3 m 2 _ ^ p 2 + m 2^/ k i + m 2 +pk + t£ 

1 n B (\/p 2 + m 2 /T\ n B (Vk 2 + m 2 /T\ 

V k 2 + m 2 \Jp 2 + m 2 ^ ' ^ 

^ f \J k 2 + m 2 a/p 2 + m 2 + fcp j ^ 



16 



V 



k l p 2 - (kpf 
(p + k) 2 



^ ^ \/fc 2 + m 2 a/p 2 + m 2 + fcpj ^ 



8 + 4 



V 



/J 



(4.66) 



rescale momenta by setting 



p = x ■ |0| and k = y 



(4.67) 



50 to polar coordinates, and introduce dimensionless variables A and A as in Eq. (|4.27j) . 
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Table 4.1: Results of numerical evaluation of the singular part of AP^^/P one _ loop for 
A = 20 and A = 100. In the limit e — > 0, the real part converges to a finite value; with 
e = 10" 8 , already four digits are stable. The imaginary part goes to zero for vanishing e. 



Finally we arrive at 



AR 



MHH 
VTT 



e 2 A 4 A" 



dx dy dt 



x 2 y 2 



8(2tt) 4 J J J Vx 2 + 4e 2 vV + 4e 2 
n B (2ttA- 3/ Vx 2 + 4e 2 ) n B (2ttA" 3/2 vV + 4e 2 

xV(t 2 - 1) 



a; 2 + y 2 + 2xyt 



4e 2 — \/ x 1 + 4e 2 \/y 2 + 4e 2 — a;yt + 



^2„,2 



I () ( x 2 + y 2 + ^(1 + t 2 ) + gVx 2 + 4e 2 v / l/ 2 + 4e 2 



(4.68) 



The condition in Eq. ()4.64j) is implemented exactly as in the case of AP™, see Eq. 
to (j4.39j) ; cos 6 has been abbreviated by t in Eq. (|4.68j) . 

As already indicated above, the term proportional to ^2^2+2^ ^ as a singularity. We 
evaluate the integral numerically by prescribing a small value for e. In Tab. 14.11 one can 
see that for decreasing e the real part converges to a finite value, while the imaginary part 
converges to zero. For our purposes, working with e = 10 -8 is sufficient for determining the 
real part; the imaginary part is ignored. The second term is regular, and the evaluation is 
straightforward by setting e = 0. The result of numerical evaluation is shown in Fig.[ 



4.3 Results 

To compare the corrections arising from the two-loop diagrams to the one-loop pressure, we 
plot the ratio of each of the contributions to AP and P on e-ioop from Eq. ()4.3|) as a function 
of the dimensionless temperature A c = 11.65 < A < 250. The results are shown in Fig. 14.61 
through Fig. 14.91 The one-loop pressure does not include the ground-state contribution 
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Eq. (|4.2jl . For the effective gauge coupling, the plateau value e = 5.1 is used for all tem- 
peratures A. Throughout most of the electric phase, this is admissible, but the logarithmic 
pole of e at the critical temperature A c is ignored. The nonlocal diagram is the dominating 
contribution for A < 100. Throughout the electric phase, the corrections arising from two- 
loop diagrams are tiny; the ratio of two-loop to one-loop contribution is 2 • 10~ 3 at most. 
With rising temperature, the two-loop contributions decrease. The contributions AP^Sjf 1 
and AP™ approach zero for large temperatures. In contrast to that, AP^t 1 becomes con- 
stant. For APyJjP we obtained only a rough estimate. For temperatures close to the phase 
transition, there is a dip in the two-loop corrections to the pressure. The minimum of the 
dominating contribution is at A ~ 27. 

The microscopic interpretation is as follows: Close to the phase transition at A c , the 
monopole mass M monopo i e oc — decreases sizeably. This increases the scattering of TLM 
modes off magnetic monopoles, with the consequence of a dip in the dominating two-loop 
contribution. With rising temperature the monopoles become massive and dilute, and scat- 
tering processes are suppressed. But as even for asymptotically high temperatures massive 
and dilute scattering centers are present, the contribution AP-f^r remains finite. 
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Figure 4.6: The ratio of APj^ and P „ c -ioop plotted for 11.65 < A < 100 and 11.65 < A < 
250. 
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Figure 4.7: The ratio of AP™ T H and P on c-ioo P plotted for 11.65 < A < 100 and 11.65 < A < 
250. 
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Figure 4.8: An upper bound for the ratio of APy^ 1 and P on e-ioop plotted for 11.65 < A < 100 
and 11.65 < A < 250. 
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Figure 4.9: The ratio of AP™ and P onc _ioo P plotted for 11.65 < A < 100 and 11.65 < A < 
250. 
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Chapter 5 

Summary and Outlook 



In this thesis, we computed the phase and the modulus of a composite, adjoint Higgs field 
relevant for the thermodynamical description of an SU(2) Yang-Mills theory. The phase 
was defined by giving a spatial and moduli space average over an adjointly transforming 
two-point function and demanding BPS-saturation. The modulus could be inferred after 
assuming the existence of an externally given scale. It was seen that the field <fi exploits 
the instantaneous long range correlations in the classical caloron solution. The nontrivial 
temporal winding was obtained only after averaging over the entire admissible part of the 
moduli space. From the BPS equation, 0's potential was uniquely deduced. 

The two-loop contribution to the pressure of SU(2) Yang-Mills theory was computed. 
There are four contributing Feynman diagrams, among them one nonlocal. This one is 
seen to be dominant. The two-loop corrections are smaller than 2 ■ 10~ 3 as compared 
to the one-loop result. Therefore the underlying picture of only very weakly interacting 
thermal quasiparticles is confirmed. The importance of the compositeness constraints in this 
process was pointed out. A microscopic interpretation in terms of TLM modes scattering 
off magnetic monopoles was given. 

We saw that the compositeness constraints rule out (at least in the large coupling regime) 
a number of diagrams which one would naively expect to contribute. It would be interesting 
to consider this for higher loop-diagrams. Possibly, from some loop order upward only a 
few classes of diagrams will survive. The effects of the two-loop corrections for the pressure 
on the evolution of the gauge coupling are currently being worked on [2H]- A review on 
the discussed approach to SU(2) and SU(3) Yang-Mills thermodynamics (containing also 
the material presented here) and especially indicating implications on particle physics and 
cosmology is available in 0. 
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Appendix A 
Appendix to Chapter 3 



A.l Notation and conventions 

In Chapter 3, the following conventions are used: SU(2) gauge field and field strength are written 
in matrix notation as 

^ = ^y and F ^ = F^^, (A.l) 
where A a are the Pauli matrices. The covariant derivative is defined as 

D M = <9 M - iA„ , (A.2) 
where the coupling constant is absorbed in the gauge field. The field strength is in matrix notation 

F„ v = A v - d v Ap - i L4 M , A v \ (A.3) 
and in components 

= dpA a v - d v Al + e abc AlAl . (AA) 



We work in a Euclidean spacetime with x = (t,x). Latin indices run from 1 to 3, Greek indices 
from 1 to 4. Upper and lower indices are not distinguished. The summation convention is implied 
if not specified otherwise. The totally antisymmetric symbol ep Upa has £1234 = e 1234 = +1. The 
four-gradient is 



8 " = aV < A ' 5 > 

derivatives with respect to time and space coordinates are written as 

® A = Jr and di = llx- 
respectively. The derivative with respect to the radial coordinate r = \x\ is written as 

d r = I . (A.7) 
The Pauli matrices are denoted by A a , they have the explicit form 

, i fo i\ x2 _ fo -A x3 ( 1 
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The Pauli matrices fulfill 



A", A = 2ie a b c X 

&ab + i£abcA 

The traces of products of Pauli-matrices are 



A a A b 
A a A b A c 



(A.9) 



S a A b + S bc X a + ie 



abc ■ 



trA a 
trA a A fe 



trA a A 6 A c 





tr A a A 6 A c A d = 2 (5 a b5 C( i + 5 a dhc - S ac ^bd) 

tr A a A b A C A d A £ = 2i (5 a b£cde + ^cd^abe ~ ^cd^abd + 5de£abc) 

The 't Hooft symbols r] and fj are defined as 



(A.10) 



Vpu 



'dpv 



+ Sap&vi 
SauSv 



5av&pi 
SauSpA 



(A.ll) 



The symbols r\ (77) are (anti) self-dual and antisymmetric in the vector indices, 
1 



jJLU 



v 

a 

fXU 

a 



'^pvafi Va/3 



(A.12) 



-v, 



up 
a 



They fulfill a number of useful relations: 



VpvVpv 


= 4<5 a6 




r IX f iV\v 


= 35fj, u 




VpvVpv 


= 12 




'HxpVxu 


= 5ab&pv 


+ ^abcVpv 


b c 
^abcVpu^KX 






b c 
e abc 7 ! pu 7 ! pX 


= 2^a 




VpuVvXVXp 


^abc ■ 





(A.13) 



MpX 



The relations ()A.13|) hold for r/ as well. Besides that, 

VpuVnX = ^pn^uX — &pX&vn + ZflUKX 



n^nnx 



OflK"l/X ^pX^UK £pvnX 

^XpuaVpa = ^pX^pu + ^pvVXp + ^PpVuk 
-ZXpvvVpa = $pXVpv + dpvVXp, + dppfjvK 



(A.14) 
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A. 2 Caloron field strength 

The SU(2) field strength in the 'nonperturbative' convention is 

F« v = d,A a u - d„A« + e^A^Al . (A.15) 
Inserting the singular gauge instanton field 

Al = -^ K d K \nU = -fj; K ^, (A.16) 
we get 

—a n(^n)-an)^n) ^ Q ijjw) - (9 K n)(^n) 

/ (fljiX&n) _ (a„n)(a A n) _ (d K n)(d M n) 
1 ^ n 2 mA n 2 n 2 
_ a (^n)(^n) _ _ a n(^^n)-2(9 K n)(^n) _ a n(^^n) - 2(d K u)(d v n) 

Vfj,u -q2 ^Iuk -q2 + ^ JT2 

For convenience, we define the quantities 



P= («p a„ d P mmzmmm, (A . W) 

ri 2 n 2 

such that the instanton field strength reads 



C = V^P ~ + • (A.19) 

The symbols (|A.18|) have the properties 

P»v = Pw and P^ = -2P. (A.20) 

The pre-potential II for a single caloron is given in Eq. ()2.1()4j) . Using the fact that it depends 
on the radial coordinate r = \x\ only (and on Euclidean time), H(r,x) = U(r,r), P and P^ v can 
be expressed in terms of the derivatives of II with respect to r and r. In particular, the Laplace 
equation 

n- 1 9 M 5 M n = o, (A.21) 

which results from the demand for self-duality, translates into 

f + f + 2 M.„. (A, 2 ) 
At an arbitrary point x = (t,x), the components of P^ are explicitly given as 

dju(x) ja 4 n(x)] 2 



Pu (x) 



n(a?) n 2 ( 



" 7 V n(x) " iPp J (A - 23) 

„ . . . d r U( X^j X^Xj ( d^H(x) d r IL(x) [d r IL(x) 

±lj [X ) = Oij r r 



rll(x) r 2 y IL(x) rTI(;r) Il 2 (2;) 
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and 

_ t^n(x)] 2 [g,n(x)] 2 
1 ] ~ n=(x) + n2(x) ' 1 j 

At x = (r, 0) , the tensor P^ (r, 0) is diagonal and its components are 
P 4i (r,0) = P i4 (r,0) =0 



a|n(r,o) ja 4 n(r,o)] 2 



P44(T ' 0) " IT(r,0) 2 n2(r,'o) (A.25) 
uv ' ; 3 J II(r,0) 

and 



A. 3 Details to section 13.2.3 



All spacetime dependent objects are to be evaluated at the same time r, so that we may, for 
simplicity of notation, suppress the time dependence and write F fiu (x) instead of F pu (r,x) and 
{0,x} instead of {(r, 0), (r, a;)} etc. Besides that, we abbreviate the time derivative with a dot 
and the radial derivative with a prime. 

In the following some expressions containing P pu and P, which we will need later, are calculated: 

nl K nl P PuM p w (x) = (v<W - <WM JwO") 

= P^ (0) P M „ (x) - P w (0) P w (as) ( A.27) 

= P^(0)P^(a;)-4P(0)P( a; ) 



^Cp^(0)^p(^) = ^4^4pA4(0)P w ( a; ) + ^ fc ^ p P nfe (0)P w ( a; ) 

^ c ii(Q) 
5 n(o) 



-P 44 (0)P ac ( a; ) - -f\l n f\ c np ^P w {x) 



-Pu(tyPac(x) - ^(e aiin + S^San) (e cnp - 5 cn 5 P 4)^|^yP At p(a;) 
-P 44 (0)P ac (x) - l^L(p ac ( x ) - S ac Pii(x) - 5 ac P u (x)) (A.28) 

-P 44 (0)Pac(x) - ~§:{Pac(x) ~ 5 ac P^(x)) 

-Pu(0)P ac (x) - ^(P 4 4(0) + 2P(0)) (P ac (x) + 25 ac P(x)) 
~(2P u (0)+P(0))P ac (x) - ^(P u (0) + 2P(0))5 ac P(x) 



, . fi(r) n /2 (r) IK . 
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P^(0) P^fx) = P 44 (0) P 44 (:z) + Pmn(0) P mn (x) 

iii(o) 



P 44 (0)P 44 (:z) 



311(0) 



Pii{x) 



/fi(o) ri 2 (o)\ /fi(r) n 2 (r)' 
^n(o) n 2 (o) i I n(r) n 2 (r) 



, in(Q) /n(r) | 2 



n /2 ( 



3n(o) n(r) n 2 (r) 



(A.30) 



^uW^{ X ) nvn*) - y u2{0) 3 n( ) ) [ 6 W( r ) 2 n(r) n 2 ( 



fj a Kp P MK (0) P w (a;) = ^ ^44(0) P 4p (*) + 77g p P mfc (0) P mp { 
= Pu(0)P4a&) ~ n a m p Pm P {x) 

= p 44 (o)p 4a (*) + i5j|p 4a ( a; ) 



_ x a ( tl'jr) tl(r)W(r) \ ( 2 ft(0) _ , ; 

r ln(r) n 2 (r) 311(0) n 2 (o)J 



ri 2 (o)\ 



(A.31) 



(A.32) 



x a xf 



r 2 ^I'np ± /ARK" J ± MP - 



n(r 



n 2 (r) i I 311(0) n 2 (o) J 



(A.33) 



Integrand of 1)3.31(1 . Inserting the result for the Wilson lines from Eq. (|3.37|) and writing the field 
strength in components, we have Eq. 1)3. 42|) . which in simplified notation reads 



trA a P M ,(0) {0,x} F^x) {x,0} 



1 

-tr 
2 



A a A b ( cos g(r) + i\ c — sing(r) ) \ d ( cos g(r) — i\ e — sin #(r) 



i 



cos^(r)-tr[A a A 6 A d ]-P^(0)Pt(^ 



- -sin 5 (r) cos 5 (r) ■ tr[A a A fe A d A e ] • - • F%,(0) F%(x) 



i 

2 

i 

2 
1 



(A.34) 



+ 1 S mg(r) cosg(r) • tr[A a A fe A c A d ] • X - ■ F%,(0) F^(x) 



+ ^ sin 2 g{r) ■tr[\ a \ b \ c \ d \ e ] 



We will calculate each of the four summands separately. Therefore we write the Lorentz-trace of 
the field components F b u (0)F d u (x) in terms of the symbols P and P VbV introduced in Eqs. (|A.18|I . 

F* v (Q)F* v (x) 

= [^ P (°) + ^« P -( ) " VLP^(0)] [n%P{K) + t p Pu P (x) - ntpP^)] (A.35) 
= 2{$ M [iV(0)iVOiO - 2P(0)P(x)] +e bdf 4 p P llK (0)P IJip (x) -ff^tpPu^P^x)] . 
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Furthermore, we need the traces of products of Pauli matrices Eqs. (|A.10|) . some of the relations 
for 't Hooft symbols from Sec. IA.ll and the expressions prepared in Eqs. (|A.27|) through (|A.33|) . 
The first of the four terms in Eq. (|A.34|) is 



icosVr-)-tr[A a A"A d ]-F^(0)^(^) 

= zcos 2 g(r).e abd F b „(0)F d Jx) 
= 2i cos 2 g(r)-t Kp P^(0)Pw(x) 

= -W,(r) 2^P) . (A.36) 

yv 1 r ^n 2 (o) 3n(o) J ^n(r) n 2 (r) J v ; 

In the third line, we used that 

e ab % K ni p P VK {$)P w {x) = e adb fj d pi/ fj b Kp P KU (0)P pp (x) = -e ab % K vi p P UK (0)P» P (x) 
vanishes. The second summand is 

- I cos g(r) smg(r) • tr[A a A b A d A e ] • ^ • F b u (0) F% 



cos g(r) sin g(r) ■ tr[A a A 6 A d A e ] • - • F b u (0) F d v {x) 
icosg{r) smg(r) ■ {5 ab 5 de - 5 ac 5 be + 5 bd 5 ae ) ■ y ■ F b u (0) F d u (x) 
i cos g(r) sin g(r) • ^ • ^35 ae [P^O) P^x) - 2P(0) P(x)] + 2e aef fj f Kp P MK (0) P pp {x) 

icos 5 (r) sin g(r) ■ ^- ■ {3[P pu (0)P pu (x)-2P(0)P(x)] - [P^(0) P^x) - 4P(0) P(x)] } 
2i cos g(r) sin g(r) ■ ^ • [p^(0)P^(x) - P(0)P(x)} 



I n 2 (o) 3n(o) J \ n 2 (r) n(r) n 2 (? 



The third summand is 



- cos 5 (r) sin 5 (r) • tr[A a A b A c A d ] • - • F^(0) F^(x) 
= i cos g(r) s\ng{r)-{5 ab 5 cd -5 ac 5 bd + 5 bc 5 ad )-^-F b v {U)F d v {x) 
= i cos g(r) sin g(r) ~'{~ M P ^(°) P A X ) ~ 2P(0)P(a:)] 

" {%n€ P + i^«(0) P w {x) + tfac^^p P-(0) P w (^)} 

= -2icosg{r) smg(r) ■ { ^P(O) P(x) + P„«(0) P w (x) 

, , x a / n 2 (o) 2 ri(o) \ f ii 2 (r) n ,2 (r) n"(r)\ lk x 

= 2icos5(r) sin 9 (r) — ^ - -ttH TT^H - 3tt^t + 2 T7^ • A -38 

yw yw r I n 2 (0 3h J \lP(r) U 2 (r) U(r) J v ; 
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The fourth and last one is 



1 . 9 / \ r, /i , A, ^ , //, on XX 

2 



J sin^r) • tr[A a A b A c A d A e ] • — • **,(0) <(x) 



zsin 2 5f(r) • (<5 a6 e c d e + (5 cd e abe - <5 ce e aM + S de e abc ) ■ ■ F b {Q) F*(x) 



r 2 [iv \ J fiv 

ism 2 g(r) ■ ( 2 £ a6c ^ ~ ^abd) ■ Ft{0) F d (x) 



ism 2 g(r) ■ j Ue abc e M ^- - 25 a ^j r) f Kp P^ K (0) P w (x) - 2e ahc ^-j]\ K rf vp P UK (0) P w (x) 



= ~^(r) ■ X ~ ■ ~ |5© - 2^P) . (A.39) 

yy ' r ln 2 (0) 311(0) j\ n(r) n 2 (r) J v 7 

Summing Eqs. (|A,36|) through (jA.39j) and using 

cos 2 a — sin 2 a = cos 2a and cos a sin a = — sin 2a (A. 40) 

yields (with the r's reinserted) 

trA a F„„(T,0) {(t,0),(t,*)} F„(T,!t) {(t, x), (t, 0)} 

= 2 . £ / tf^ 2 fi^oA f / ri(r, r) n-(r, r) fr(r, r) \ 



n 2 (r,0) 3H(t,0) I [ "1 n 2 (r,r) n(r,r) I 

+ sin(2a(r r)) f 2 5^ - 2^1 + ^) _ 1 (A 41) 

+ S 1 51 ' jj { n 2 (r,r) n 2 (r,r) + H(r,r) H(r,r) J f " 1 ' 



Appendix B 
Appendix to Chapter 4 



B.l Details to section 14.2 



In this appendix, the contraction of Lorentz and color indices appearing in the expressions for the 
Feynman diagrams Figs, 14.21 14.41 l4~o1 is performed. For diagrams containing TLM fluctuations and 
hence P^ v as defined in Eq. (|4.11jl . the following formulae will be useful: 



P»p^=pfc-M^ (B.l) 



Concerning AP HH . The color indices a, b, c and d take the values 1,2 only because we chose the 
particles propagating in the loops to be TLH-modes. Index / is summed over 1,2,3. 



[e adf e }hc (g^g pcT - g» p g va ) + e abf e fdc (g^g vp - g^g™) + e acf e fdb {g^ g vp - g^g° p )] 

= e acf e acf - sTfT - 9^9 vp + iTsT) • (<V " ^) ■ (W - ^) (B.2) 

= 4 fl2-34-34 + ^-^ 

\ m z m z mr J 



Concerning AP . Since k and p are associated with a TLM-mode and a TLH-mode respectively, 
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we have to set a = b = 3 and sum over c, d = 1, 2. Eqs. (|B.1|) have to be used. 

[e adf e fbc (g pu g p ° - g^g va ) + e abf e fdc (g p ° g» p - gTsT) + e acf e fdb (g^f - sTf) 

■(Sod) (9pa-^)-(-Sa,)P^(k) 



eacfSacf iaTgr - w - <r<r + <rf p ) • [g P , 

Hg^g^-g^gn-igp*-^ 



PpPu 

m 2 



PL(k) 



(B.3) 



, , „ , P 2 (Pfe) 2 

-6 + 2^ H ^ TT2 



Concerning AP MHH . The momentum g is associated with a TLM fluctuation, and the TLH 
fluctuations carry momenta p and k. Thus we have f,g = 3 and a,b,c,d= 1,2, and hence 



2 2 



£acf£bdg$ab$cd,dfg ~ ^ ^ £ac3£ac3 — 2 . 



(B.4) 



a=l c=l 

The Lorentz structure resulting from the propagators and polarization tensors is 

g w (k - P ) x + g px (q - kY + g x »(p - q) p ] \g va (p - k) K + g aK {k - q) u + g KU (g - pf] 



■ [g^u 



PjlPv 

m 2 



g P a 



kpk a 



pIM- 



(B.5) 



Contracting the Lorentz indices and inserting P\L(q) by using the rules (jB.l|) yields 



2 k 



(kg)' 



+ 2k 2 - Akq - 4pq + 2p 2 + 4q 2 



+ 



P 



)- 2 (/-^) +6 ( fe p-«) 

2 ^)(p 2 + q 2 )-,(k P -^)kp 



(B.6) 



4(A;p) 2 + 4(kp)(kq) - 2{kq) 2 + 4(fcp)(pg) - 2{pqf 



+ 



1 



(M 5 



(kqf 



{pqY + 2 [kp 



(kq)(pq) 



(kq)(pq) 



Now we may use energy momentum conservation and set q = —p — k. Moreover, we notice that 

,2 , i„„\2 



k 2_K±kpY =p 2 



(k+p) 2 * (k+p) 2 

This leads to the following: 
{kp) 2 - k 2 p 2 



p 2 + kp) 2 (k 2 + kp)(p 2 + kp) , k 2 p 2 - (kp) 2 

fcp — 



(fc + p) 



(fc + p) 5 



10- 



(k + p) 



+ 10fc 2 + 16A;p + 10p 2 



(B.7) 



(B.8) 



1 



2^2 



+ ( -2k 4 - 8k 2 kp - I6(kp) 2 - 8kpp 2 - 2p 4 - %^-(Zk 2 + 8kp + 3p 2 ; 



(k + p) 



1 (fep)^ k 2 p 2 ^ 4 + Ak2kp + + ^ + 4fc ^ 2 + ^ 



+ 



m 4 (fc+p) 2 
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k 



fj,, a = 3 




u, b = 3 



Figure B.l: The tadpole contribution to the TLM mode polarization tensor. 



We are interested in only one contribution from this diagram, namely APyJ^f 1 . Fortunately, its 
calculation includes the two on mass shell conditions p 2 = m 2 and k 2 = m 2 . They reduce the 
former expression to 

M (rf.jf).*VlM f 8 + 4<^V (B.9) 



m? J (p + k) 2 



B.2 Electric screening mass for TLM modes 

There are two nonvanishing contributions to the polarization tensor Jl fJiU (p) for TLM modes on one- 
loop level. They are depicted in Figs . IB . H and lB~2l We are especially interested in the component 
IToo for vanishing external momentum, more precisely Hqo(pq = 0, \p\ — > 0). 



B.2.1 Tadpole contribution 

According to the Feynman rules (see Sec. I4.1|) . the tadpole diagram Fig. IB.ll corresponds to the 



expression 



iC P >) = - 



d 4 k 

(2tt)0 



kpkfj 



k 2 



+ 27r5(k 2 -m 2 )n B (\k \/T) 



■ (-ie 2 ) ■ [e eab e ecd (g^g™ - gT <f p ) + e eac e edb {sTgT - g^g pa ) 
+e ead e ebc (g^g ap -g^g a n} 



(B.10) 



d 4 k 

W) A 



k 2 



m- 



+ 



k^k u 



k 2 



+ 2n5{k 2 -m 2 )n B (\k \/T) 



When integrating the loop momentum k, the following constraints have to be taken into account: 
1. For both thermal and vacuum contribution the constraint on the center-of-mass energy in 
vertices, 



\(p + k) 2 \ < 



(B.ll) 



2. For the vacuum contribution additionally the constraint on the offshellness of quantum 
fluctuations, 



\k 



m 2 \ < 



(B.12) 



As for the calculations of two-loop diagrams in Sec. 4, massive vacuum fluctuations for e > \ are 
forbidden by Eq. (jB.12|) . i. e. the vacuum contribution vanishes for e > ^. Moreover, for e > \ and 
external momentum p = 0, the thermal part also vanishes because in this case 



Up + kf 



m 2 = (2e) 2 H 2 > 



(B.13) 



and Eq. (|B,11|) cannot be satisfied. Thus we have no contribution from the tadpole for e > 1/2. 
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fi, a = 3 




b = 3 



Figure B.2: The nonlocal contribution to the TLM mode polarization tensor. 



B.2.2 Nonlocal contribution 



The nonlocal contribution to the polarization tensor is depicted in Fig. IB.2l In formula, it reads 



n nonl» = ^ 



1 f d 4 k 



2i J (2vr) 4 



e • s„ 



g^(-p - k) x + g px {k -p+ky + g Xp (p - k + pf 



■ £dbf 



g av {-k -p) k + g UK {p + P- k)° + g Ka {-p + k + k) u 



(Scd) ( 9pa 



kpk a 



k 2 



+ 2it5{k 2 -m 2 )n B (\k \/T) 



■ (S e f) 9\ 



(P- k) x (p- k) K 



(p — k) 2 — m 2 
d A k { 



te 



+ 2vr 5{(p - k) - m 2 ) n B (bo - fo| /T) 



2^-2^-2^ + 4^ * {kp) 
m 2 m 2 



i^- + 5^ + 2^-^],r 



+ (10 + 2-^-2-^-34 + 41 ^ + f- 2 " 3 4 + ^T- h — ) 

m z m z m z rrr I \ m z m z 



(kp) 2 i p 2 



+ 



m z m z 



+ 2vr5(fe 2 -m 2 )n b (| A: |/r) 

+ 2vr 5{(p - kf - m 2 ) n B (bo - fco| /T) 



k 2 — m 2 
i 



(p — k) 2 — m 



(B.14) 



Again, all contributions containing a vacuum TLH fluctuation vanish for e > 1/2. The only 
nonvanishing contribution is the case of both fluctuations being thermal; it reads 



n 



nonl. 
therm. 



(p) = ie 2 



d 4 k 



{kpf 



•2k p - l^-r- + 7p 2 - 4 ) sT + ( 12 - 2^- - 3^ + ^ ) k»k> 



kp 
m 2 



1 4 
P p 



mr 



+ -6 + 4 



kp p 2 kp 



mr 



)( fc v + fc v) + (-5+^! + 4)^} 



S(k 2 - m 2 ) n B (\k \ /T) ■ 5((p - k) 2 - m 2 ) n B (|p - *o| /T) . 



(B.15) 
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We now set po = (keeping |p| / 0) and get 



lC nl >o = 0,p) 



therm 

= ie 2 



d 4 k 
(2vr) 2 



(kp) 

-2fcp-4^f 

rn * 



+ -6-4 



kp p 2 kp 

m 2 m 4 



m^ 1 m 



^2 ) pV 



5(fe 2 - to 2 ) • 5((p - A;) 2 - to 2 ) • [n B (|fc | /T) f 



The two 5-functions can for p ^ be rewritten as 



5(A; 2 - m 2 ) • 5{{p - kf - to 2 ) = 5(k 2 - m 2 ) ■ 5(p 2 - 2pk) 



(B.16) 



1 



6 (k - \Jk 2 + to 2 ) • (5 (p 2 , - p 2 - 2p V^ + ro 2 + 2pfc 



2\/fc 2 + m 2 

+ 6 (k + v / fc 2 ~Tm 2 ) • 5 (p 2 , ~ P 2 + ZpoVk 2 + m 2 + 2pfc^ 

- |p| 2 -2poV|fc 2 + ^ 2 + 2|p||fc| 



2vW + 



TO' 



|W + #F+^ -<MPo- |p| 2 + 2poV\ k \ 2 + m 2 + 2|p||fc| cos6> 



(B.17) 



In the case po = (with |p| ^ 0) the former expression reduces to 



5{k 2 — m 2 ) • 5{{jp — k) 2 — to 2 ) 



2^ 



TO' 



(5 (& - y/\k\ 2 + to 2 ) • 5 (2|p||fe| cos# - |p| 2 ) 



(B.18) 



+ <5 (fc + ^|fc| 2 + TO 2 ) • 5 (2\p\\k\ cosO- \p\ 2 ) 
1 -/ .. \P\ 



• 5 cos i 



4\p\\k\^/\k\ 2 + m 2 " V"' 2 |fc| 
Inserting Eq. ()B.18|) into Eq. (|B.16|) yields 



5 fc - V|fc| 2 + TO 2 )+5 (k + ^\k\ 2 + m 2 



ni nl .(Po = o )P ) 

therm. 

2 /-co 



(27T) 



2 i iMivQ 

oo ./-p/2 



A; 2 ^ / dU 



n B 



{Vk 2 + m 2 /T) 



| (-8p 2 - + f 12 + 4-^ + k^k v 

I \ to y \ to to y 



(B.19) 



+ -6-2 



p 1 p 
m 2 2 to 4 



<5(cos I 



2k' 



ApkV k 2 + to 2 



5 ( A;q — ytf + m?) + 5 (ko + \Jk 2 + rri 
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Only the terms proportional to g^ v or k^k u contribute to II . The result is 

n° n ° on i. (po = o, p ) 

therm. 



!!L(ll + 4P + P) r dk ky / k 2 + m 2 n 2 B ( Vk 2 + m 2 /T) . 
47r V p m 2 m 4 J J p/2 V ) 



In the limit |p| — > this expression diverges, 

3ie 2 1 



n°° onl (po = 0,p) ► — A r dkk Vk 2 + m 2 n% (Vk 2 + m 2 /T) . (B.21) 

therm. IpHO VT |p| J Q \ J 
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